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ABSTRACT 

We introduce a novel formalism to investigate the role of the spin angular momentum of 
astrophysical black holes in influencing the behaviour of low angular momentum general 
relativistic accretion. We propose a metric independent analysis of axisymmetric general rel- 
ativistic flow, and consequently formulate the space and time dependent equations describing 
the general relativistic hydrodynamic accretion flow in the Kerr metric. The associated sta- 
tionary critical solutions for such flow equations are provided, as well as the stability of the 
stationary transonic configuration is examined using a novel linear perturbation technique. 
We examine the properties of infalling material for both the prograde as well as the retrograde 
accretion as a function of the Kerr parameter at the extreme close proximity of the event hori- 
zon. Our formalism can be used to identify a new spectral signature of black hole spin, and has 
the potential of performing the black hole shadow imaging corresponding to the low angular 
momentum accretion flow. 
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1 INTRODUCTION 

Astrophysical black holes are the terminal states of the gravitational collapse of massive celestial objects that can be conceived as singularities 
in space time censored by a mathematically defined 'one way barrier' - the event horizon, and are not amenable to any direct physical 
observation. As a result, their presence can only be realized through the gravitational influence they exerts on the matter accreting onto 
those objects. The infalling matter inevitably plunges through the event horizon on a relativistic scale of velocity. Given a set of physically 
reaUzable outer boundary conditions, such accretion eventually manifests transonic properties to obey the aforementioned inner boundary 
conditions (Liang & Thomson 1980). Subsonic at a large distance, accretion thus reaches the event horizon supersonically. 

The hypothesis that most (if not all) of the supermassive black holes and the stellar mas black holes powering the active galactic 
nuclei and the galactic micro quasars respectively possesses non zero value of spin angular momentimi has gained widespread currency in 
recent times (Miller, Reynolds, Fabian, Miniutti & Gallo 2009; Kato, Miyoshi, Takahashi, Negoro &. Matsumoto 2010; Ziolkowski 2010; 
Tchekhovskoy,Narayan & McKinney 2010; Daly 2011; Buliga, Globina, Gnedin, Natsvlishvili, Pitrovich, & Shakht 2011; Reynolds et al. 
2011; McClintock et al. 2011; Martinez-Sansigre & Rawlings 2011; Dauser, Wilms, Reynolds & Brermeman 2010; Nixon, Cossing, King & 
Pringle 201 1 ; Tchekhovskoy & McKinney 2012). The black hole spin plays a deterministic role in influencing various characteristic features 
of the the dynamical and spectral features of accretion and related phenomena in the characteristic black hole metric - the energy extraction 
from a Kerr black hole through the Blandford Znajek Mechanism (Garofalo 2009a; Das & Czemy 2012a; Tchekhovskoy & McKiimey 
2012), the spin dependence of the black hole shadow imaging (Falcke, Melia & Agol 2000; Takahashi 2004; Huang, Cai, Shen & Yuan 
2007; Hioki & Maeda 2009; Zakharov, Paolis, Ingrosso & Nucita 2012; Straub, Vincent, Abramowicz, Gourgoulhon, & Paumard 2012), 
various evolutionary properties of the normal and the active galaxies (Garofalo 2009b; Garofalo, Evans, & Sambuma 2010; McNamara, 
Rohanizadegan & Nulsen 2011), and the QPO phenomena associated with the extra-galactic sources (Czemy, Lachowicz, Dovciak, Karas, 
Pechacek, & Das 2010; Das & Czemy 201 1), to mention a few. 

The aforementioned issues necessitate investigation of how the black hole spin angular momentum, i.e., the Kerr parameter a, influences 
the dynamical and the radiative behaviour of the general relativistic transonic accretion at the close vicinity of the event horizon of a rotating 
black hole. The central objective of our work presented in this paper is to investigate what properties of the accretion flow are the principal 
attributes of the black hole spin at the extreme close proximity of the event horizon of a rotating black hole. 

To accomplish such task, we conduct a detailed and multi-step investigation of the transonic properties of general relativistic axisymmet- 
ric hydrodynamic inviscid accretion of low angular momentum as realized on the equatorial plane of the Kerr metric using the Boyer Linsquist 
(Boyer & Lindquist 1967) co-ordinate. We begin with a general prescription where we consider a four dimensional stationary axisymmetric 
manifold with two commuting Killing vector fields and subsequently construct the general relativistic Euler and the continuity equation from 
the appropriate energy momentum tensor. Quite interestingly, we have been able to demonstrate, using certain symmetry arguments, that for 
the three dimensional submanifold the fluid equations are separable using analytical scheme and the corresponding flow velocity components 
can completely be determined once the equation of state is specified. We thus formulate a general framework for studying the equations for 
fluid flow in rotating black hole spacetime. 

The transonic flow properties in phase portrait, however, can not be determined completely analytically for certain issues which will be 
elaborated in the subsequent sections. For such purpose, we demonstrate the equivalence between the set of differential equations describing 
the space gradient of the stationary flow velocity and that of the acoustic perturbation (stationary barotropic sound speed), and the first order 
autonomous dynamical systems. The emergence of the multi-transonic behaviour manifests through the critical point analysis. Such multi- 
transonic accretion solution, as we will see in the subsequent sections, may contain stationary shock, properties of which are obtained by the 
explicit solution of the general relativistic Rankine-Hugoniot conditions. The properties of the post shock flow are then studied as a function 
of the black hole spin - the Kerr parameter a. The post shock flow solution are then followed upto a sufficiently close proximity of the event 
horizon to demonstrate how the terminal values of the shocked accretion variables are influenced by the black ole spin angular momentum, 
and the consequences of such dependence are discussed in detail. 

The entire procedure to the study the spin dependence of the behaviour of accreting matter close to the event horizon as described above 

is based on the stationary solutions of the differential equations describing the accretion phenomena. Along with the understanding of the 
transonic behaviour of the stationary flow solutions, it is equally important (rather necessary) to ensure the stability of such stationary con- 
figurations. Criteria for such sustained stability can be examined by studying the time evolution of a linear acoustic like perturbation (applied 
around the stationary configuration) in the full time dependent flow equations. The existence of the stable stationary transonic solution is 
associated with the non-divergent linear perturbation of the aforementioned kind. In this work, we develop a novel linear perturbation scheme 
applicable to the axisymmetric potential flow as realized on the equatorial plane of the Kerr metric. We perturb the velocity potential corre- 
sponding to the advective velocity of the low angular momentum accretion considered in our work and demonstrate that such perturbation 
does not diverge in any significant astrophysically relevant time scale. We thus formally establish the consistency of the formalism, which, 
for the first time in the literature, has been introduced in the present work to study the black hole spin dependence of the terminal behaviour 
of accreting matter at the extreme close vicinity of the event horizon of a Kerr black hole. 
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2 MULTI-TRANSONICITY IN BLACK HOLE ACCRETION: RETROSPECTIVE AND CONTEMPORARY ASPECTS 

For accretion onto astrophysical black holes, the transonicity is characterised by a transition from the subsonic state (M < 1, where M is 
the Mach number of the flow) to the supersonic state (M > 1), or vice versa. For the present work, the Mach number M is considered to be 
the local radial Mach number for stationary transonic accretion solution. For our purpose, we define A4 to be the ratio of the local advective 
velocity u (defined in the subsequent section) and the local speed of the propagation of the acoustic perturbation (local barotropic sound speed) 
Cs as defined in the subsequent section. Such transition may be a regular one through the sonic point and is associated with the transition 
of M < 1 — > M > 1 type or may be a discontineous one through a stationary shock and is associated with the M > 1 — > M < 1 
type transition. As a consequence of the fact that the non linear equations describing the steady, inviscid stationary axisymmetric flow can 
be tailored to form a first order autonomous dynamical system, the physical transonic accretion solution for the stationary axisynmietric flow 
can formaly be realized as critical solution on the phase portrait spanned by M and the radial distance measured from the event horizon - 
see, e.g., Chaudhury, Ray & Das (2006); Goswami, Khan, Ray & Das (2007) and references therein. 

For low angular momentum sub-Keplerian accretion, such transonic features may be exhibited more than once on the phase portrait 
of the stationary solutions. Such multi-transonicity as well as the resulting shock formation phenomena for axisymmetric accretion un- 
der the influence of various post Newtonian potentials, mainly, under the influence of the Paczyhski and Wiita (Paczyhski & Wiita 1980) 
pseudo-Schwarzschild potential ^, has been widely studied in the literature, see e.g., Liang & Thomson (1980); Abramowicz & Zurek 
(1981); Muchotrzeb & Paczyhski (1982); Muchotrzeb (1983); Muchotrzeb-Czerny (1986); Chakrabarti (1989); Abramowicz & Kato (1989); 
Abramowicz & Chakrabarti (1990); Das (2002); Das et al. (2003); Fukue (2004) and references therein. 

Since a regular stationary accretion solution can not encounter more than one transonic point, multi-transonicity implies a particular 
flow configuration with three critical points where two transonic solutions through two different saddle type critical points are connected by 
a discontineous staionary shock transition. The inner boundary condition imposed by teh event horizon indicates that such a combined multi- 
transonic shocked solution originates from a large distance as a subsonic flow and encoimters the outermost saddle type sonic point to become 
supersonic for the fist time. Subjected to the appropriate initial boundary conditions, such supersonic flow makes a M > 1 — > M < 1 type 
discontinous transition through a staionary shock and the shock induced subsonic flow becomes supersonic again at the innermost saddle 
type sonic point. 

One expects that a shock formation in black-hole accretion discs might be a general phenomenon because shock waves in rotating 
astrophysical flows potentially provide an important and efficient mechanism for conversion of a significant amount of the gravitational energy 
into radiation by randomizing the directed infall motion of the accreting fluid. Hence, the shocks play an important role in governing the 
overall dynamical and radiative processes taking place in astrophysical fluids accreting onto black holes. The hot and dense post shock flow 
is considered to be a powerful diagnostic tool in understanding various astrophysical phenomena like the spectral properties of the galactic 
black hole candidates ((Chakrabarti & Titarchuk 1995)) and that of the supermassive black hole at our Galactic centre ((Moscibrodzka, 
Das & Czerny 2006)), the formation and dynamics of accretion powered galactic and extra-galactic outflows ((Das & Chakrabarti 1999; 
Chattopadhyay & Das 2007; Das & Chattopadhyay 2007)), shock induced nucleosynthesis in the black hole accretion disc and the metalicity 
of the intergalactic matter ((Mukhopadhyay 1999) and references therein), and the origin of the quasi-periodic oscillation in galactic and 
extra-galactic sources ((Sponholz & Molteni 1994; Das, Rao & Vadawale 2003; Okuda et al. 2004, 2007; Czemy, Lachowicz, Dovciak, 
Karas, Pechacek, & Das 2010; Das & Czerny 2011) and references therein), to mention a few. 

The idea of shock formation in black hole accretion flow has however contested by some authors (see, e,g,. Narayan, Mahadevan & 
Quataert (1998) and references therein for a review). Nevertheless, the issue of not finding shocks in such works perhaps lies in the fact that 
only one sonic point close to the black hole may usually be explored using the framework of the shock free advection dominated accretion 
flow solutions. 

At this point we would like to emphasize that the concept of low angular momentum (allowing the multi-transonicity and the shock 
formation) is not a theoretical abstraction. Rather for realistic astrophysical systems, such sub-Keplerian weakly rotating flows are exhibited 
in various physical situations, such as detached binary systems fed by accretion from OB stellar winds ((lUarionov & Sunyaev 1975; Liang 
& Nolan 1984)), semi-detached low-mass non-magnetic binaries ((Bisikalo et al. 1998)), and super-massive black holes fed by accretion 
from slowly rotating central stellar clusters ((lUarionov 1988; Ho 1999) and references therein). Even for a standard Keplerian accretion disc, 
turbulence may produce such low angular momentum flow (see, e.g., (Igumenshchev & Abramowicz 1999), and references therein). 

Multi-transonicity in black hole accretion has been addressed using the general relativistic framework as well. The legacy of the pioneer- 
ing contributions by Bardeen, Press & Teukolsky (1972) and Novikov & Thome (1973) to study the general relativistic axisymmetric black 
hole accretion in the Kerr metric followed two different avenues, quite often in a non overlapping fashion. One school of thought essentially 
studied the transonic accretion without paying much attention to the appearance of the multi-transonicity and the formation of shock, but 
rather putting emphasis on other crucial behaviours of the flow, see, e.g., Lemos & Letelier (1994); Riffert & Herold (1995); Abramowicz, 
Chen, Granath, & Lasota (1996); Pariev (1996); Peitz & Appl (1997); Gammie & Popham (1998); Popham & Gammie (1998); Takahashi 
(2007); Sadowski (2009) and references therein. 

The appearance of the multiple critical points nn general relativistic flow onto a spinning black hole was observed and consequently 

1 This one is believed to be the most effective black hole potential introduced in the literature tUl date, see, e.g., Artemova et al. (1996); Das (2002) and 
references therein, for further detail. 
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the formation of the standing shock has been conjectured in the alternative set of (sometime contesting the aforementioned category of work 
dealing with accretion flow without the appearance of shock transition) literature, with the main motivation to explain the spectral state by 
incorporating the physics of the post shock accretion flow, as already mentioned in the previous sections. The profound work by Fukue (1987) 
is credited as the first ever consistent attempt to provide the complete formalism for the shock formation in a multiply critical accretion flow 
in the Kerr metric, although it is worth mentioning that even before Fukue (1987), multiplicity in the critical points for the general relativistic 
axisymmetric flow was addressed (Lu 1985, 1986) without mentioning the issue of the shock formation. By revisting the concept of the 
Keplerian circular motion for rotating fluids in general relativity Lu, Yu & Young (1995) intuitively explained certain issues related to the 
shock formation for multi-transonic accretion onto a Kerr black hole. 

The full general relativistic formalism introduced by Lu (1985, 1986) and Fukue (1987) was followed by Chakrabarti (1996a,b) where 
a non relativistic calculation for the shock formation for accretion (in the hydrostatic equilibrium in the vertical direction - the flow model 
we consider in our persent work) and other issues were incorporated within the relativistic framework and some of the results valid for the 
isothermal flow had directly been applied to study the polytropic flow without adequate justification. Some of such discrepancies has recently 
been modified by Mondal (2010). Chakrabarti (1990); Lu et al. (1997a,b); Lu & Yuan (1998); Lu & Gu (2004) used the general relativistic 
shock condition to study the multi-transonic flow for the conical model ^. Peitz & Appl (1997) studied the general relativistic accretion 
for multi-transonic flow but the shock issues were not dealt in detail. While all the above works dealt with the polytropic accretion, shock 
transition in general relativistic isothermal flow was discussed in Kafatos & Yang (1994); Yang & Kafatos (1995); Yuan, Dong, & Lu (1996). 
Shocked accretion for MHD flows in Kerr geometry has also been studied by various authors (Takahashi, Rilett, Fukumura, & Tsuruta 2002; 
Takahashi, Goto, Fukumura, Rilett & Tsuruta 2006; Fukumura, Takahashi, & Tsuruta 2007). 

Meanwhile, it was realized that it is instructive to incorporate an expression for the flow thickness for flow in hydrostatic equilibrium in 
the vertical direction such that the corresponding flow equation will remain non singular on the horizon. Both the thin accretion disc as well 
as the quasi-spherical flow structure can be accommodated using such a disc height. Abramowicz, Lanza & Pervival (1997) provided such 
an expression for the general flow structure. The disc height introduced by Abramowicz, Lanza & Pervival (1997) had further been modified 
to study the multi-transonic flow structure around Kerr black holes in Barai, Das, & Wiita (2004); Goswami, Khan, Ray & Das (2007); Das 
& Czerny (2012b). 

Owing to the strong curvature of space time close to the black hole, accreting fluid is expected to manifest extreme behaviour just before 
plunging into the event horizon. The spectral signature of this tremendously hot ultra fast matter with its characteristic density and pressure 
profile is expected to provide the key features of the strong gravity space time to the close proximity of the event horizon. A detail study of 
the role of the black hole spin angular momentum in influencing the dynamical features of the transonic matter close to the event horizon is 
thus very important to perform to understand the salient features of the general relativistic black hole space time, and, in turn, to understand 
the physical properties of the Kerr metric itself. Gammie & Popham (1998) and Popham & Gammie (1998) were the first to make attempt 
to understand the flow properties close to the black hole by studying the general relativistic optically thin advection dominated accretion 
flow (ADAF) in the Kerr metric. Later on, Becker & Le (2003) applied the method of post Newtonian asymptotic analysis to investigate the 
properties of the inner region of ADAF to obtain their results that has been argued to be in agreement with the relativistic flow description. 
Subsequently, Barai, Das, & Wiita (2004) studied the influence of black hole spin in determining the properties of the accretion variables 
sufficiently close to the event horizon for multi-transonic flow, although the shock conditions were not taken into account in their work. It 
has recently been demonstrated that the multi-transonicity can only be realized through the presence of a standing shock since a smooth flow 
can never make more than one regular sonic transition Das & Czerny (2012b). 

Motivated by the aforementioned issues, in this work we would like to study the behaviour of the low angular momentum multi- 
transonic shocked accretion extremely close to the black hole event horizon. The present work differs from all other previous works on 
general relativistic accretion, including Gammie & Popham (1998); Popham & Gammie (1998) and Barai, Das, & Wiita (2004) since for 
the first time in literature, not only a multi-transonic shocked flow has been studied at the close viscinity of the event horizon to understand 
the role of the black hole spin angular momentum on determining the salient features of such flow, but also a complete description of the 
linear perturbation analysis has also been provided which ensures the stablity of such accretion solutions. In addition, a formal analytical 
description for the general fluid flow in axisymmetric black hole space time has also been provided in our work. 

As an obvious extension of Barai, Das, & Wiita (2004). we introduce the standing shock by solving the relativistic Rankine-Hugoniot 
conditions, and study the behaviour of the post shock flow to the close vicinity of the horizon. We then compare such results with the 
hypothetical flow solutions for which the flow would not pass through a shock (and hence would behave like a mono-transonic flow passing 
through the saddle type outermost sonic point formed at a large distance from the black hole event horizon) for the same set of initial 
boundary conditions describing the flow. This allows us to understand whether the shock formation phenomena can alter the dynamical and 
thermodjfnamic state of matter extremely close to the event horizon and whether such change may show up through the spectral properties 
of the black hole candidates. 

From recent theoretical and observational findings, the relevance of the counter rotating accretion in black hole astrophysics is being 
increasingly evident (Dauser, Wilms, Reynolds & Brenneman 2010; Nixon, Cossing, King & Pringle 2011; Tchekhovskoy & McKinney 



^ The conical model for the accretion was first introduced in Abramowicz & Zurek (1981) for flow under the influence of the Paczydski & Wiita (1980) black 
hole potential. 
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2012). It is thus instructive to study whether the characteristic features of the terminal values of the accretion variable for the prograde flow 
differs considerably form that of the retrograde flow. To the best of our knowledge, for the first time in literature, our work presents the 
detailed spin dependence of the terminal behaviour of infalling matter for retrograde accretion onto a Kerr black hole using the complete 
general relativistic framework. 

We, however, do not explicitly consider the viscous transport of the angular momentum and the specific angular momentum of the 
accretion flow has been taken to be invariant. Reasonably large radial advective velocity for the slowly rotating sub-Keplerian flow implies 
that the infall time scale is considerably small compared to the viscous time scale for the flow profile considered in this work. Large radial 
velocities even at larger distances are due to the fact that the angular momentum content of the accreting fluid is relatively low (Beloborodov 
& Illarionov 1991; Igumenshchev & Beloborodov 1997; Proga & Begelman 2003). The assumption of inviscid flow for the accretion profile 
under consideration is thus justified from an astrophysical point of view. Such inviscid configuration has also been address by other authors 
using detail numerical simulation works (Proga & Begelman 2003; Janiuk, Sznajder, Moscibrodzka, & Proga 2009). 



3 METRIC INDEPENDENT FORMULATION FOR GENERALIZED FLOW VELOCITY FOR AXISYMMETRIC 
ACCRETION 

As already mentioned, low angular momentum non self-gravitating axisymmetric inviscid accretion will be considered in our work. To 
begin with, consider a generic (3+1) stationary axisymmetric space-time endowed with two commuting Killing fields, within which the the 
back reactionless three dimensional general relativistic fluid equations will be looked upon. In such a space-time, the combined Euler and 
continuity eqauations take the form 

^Mv^- + £iv^p(5'"'+^;V)=0, (1) 
P 

where is the time like fibre bundle (a tangent vector field in the present context) defined on the manifold constructed by the family 
of streamlines. The normalisation condition corresponding to the velocity vector field is taken to be v^v^^ = —1. Cs is the speed of 
propagation of the acoustic perturbation embeded in the accreting fluid and p is the local rest mass energy density of the fluid. For a single 
temperature fluid, p can be replaced by the particle number density. 

The space-time considered is a stationary axisymmetric manifold of dimension four endowed with two killing vector fields ^'^ and 
so that 

V(,e.)=0 = V(^^„). (2) 

The locally timelike killing field f ^ (of norm ?) is the generator of staionarity whereas the locally spacelike killing field 0'' (of norm ip) with 
closed spacelike integral curves generate the axisymmetry. It is usually not possible to obtain any orthogonal basis for the space-time of our 
consideration since iy.4>'^ for statioirry axisymmetric space- time. We would intend to specify an orthogonal basis using which the space 
time metric can directly be expressed. To accomplish such task we first define 

:= Cm - ll^-^M = i^- ^<^M, (3) 
so that generically T^t/)'^ = 0. Norm of can thus be expressed as, 

T^T^ = -(-<;2+,2^2)=_^2^ (4) 

where is timelike and vj^ is positive. The metric element can now be expressed in the orthogonal bases as follows 

OtJii^ = -■cj'^T^T^ + (p~^(j)^(j)i, + R~^Rfj,R„ + 1?~^6^6„, (5) 

{R'', 0'^} being the spacelike basis vectors orthogonal to {T'^, 0'^}. For a stationary axisymmetric space-time the hypersurface ?^ = 

defines an ergosphere rather than the horizon, c,'^ is negativew inside the ergosphere since is spacelike in that region. On the otherhand, a 
compact zu^ = hypersurface defines a killing horizon (e.g. a blackhole event horizon in our case). This can be demonstrated by constructing 
the null geodesic congruence on such a surface. 

The aforementioned formalism is valid for a very general kind of stationary axisymmetric space-time, which includes, but certainly not 
limited to, the space-time defined by the Kerr family of solutions. With reference to eq. (5), the normalisation condition for velocity vector 
field may be expressed as 

f " -2 2.-2 2, D-2 2 , 0-2 2 , 

V V Qi^u = -f Vq +(fi Vi+ R V2 +V V3 = -1, (6) 

where wo = V/j.T'^ etc. are scalars. Contracting the equation (1) with C we obtain 

v'^V^iv'^^.) + ^ [C^VpP + {i^v^)v-V.p] = 0, (7) 
where v'^v'^V = ^v^v^V^^^^-, = is ensured by virtue of the killing equation, i.e. eq. (2). Through similar procedure we also obtain 

v^V^wi + -j [0''V^p + viv^V^p] = 0. (8) 
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Note that all the differential terms appearing in Eqs. (7-8) involve partial derivatives only since v^^^ and v\ are scalars. 

Since we consider the staionary, axisymmetric flow in three dimensions, all the directional partial derivatives with respect to cj)'^ 
vanish to yield, 

from eq. (8), R being a parameter along _R^. Integration of eq. (9) provides 



VI =^exp [- j -j'^Rj, (10) 

A being a constant to be evaluated from the initial boundary conditions. 

In a similar fashion, eq. (7) provides the expression for w^^'^, which will formally be same as v\ upto an integration constant, since 
is a killing field. One thus finds, 

VO = V^^'^ - tVp,4>^. (II) 

Substitution of vq from eq. (11) and vi from eq. (10) into eq. (6) provides the expression of v-z- 

In this section we thus provide a general formalism for evaluating all ralavant bulk velocity components of a roatating accretion flow in 
a most general axisymmetric space-time, {vo, vi, V2} are, however, the general solutions and evaluation of their specific numeriacal values 
for a particular flow configuration in a predetermined (by the accreting black hole) metric is a rather involved procedure since the integration 
constant appearing in the expressions of {vo, vi,V2} can be evaluated only if the appropriate set of initial boundary conditions are provided. 
It is also requiered that sound speed as well as the rest mass energy density is to be known a priori to find the specific values of {vi}. For our 
particular purpose, the axisymmetric space-time metric is of Kerr type and {vi} = {vt,Vr, ve, v^}, for which, the initial boundary conditions 
cannot be evaluated analytically for barotropic equation of state and for a certain geomatric configuaration of the accreting fluid. Cs and p are 
also unknown a priori. 

Our aforementioned formalism for finding the general solution of {vi} is thus useful for the flow cinfiguration with known value of 
{cs,p} and initial boundary conditions. Axisynmietric low angular momentimi accretion onto an astrophysical black body however costitutes 
a complex gravitational system for which such predetermined informations are not available in general. In subsequent sections, we thus plan 
to develop a metric specific formaUsm to imderstand the spatial velocity profie of the stationary axisymmetric flow. 



4 SPACE-TIME METRIC AND THE CONSERVATION EAQUATIONS 

Hereafter, any relevant radial distance will be scaled in units of GMbh /(? and any relevant velocity will be scaled by c, where G, Mbh 
and c are universal gravitational constant, mass of the black hole and speed of light in vacuum respectively. G = c = ^bh = lis adopted. 
The allowed domain of a, the Kerr paramter is taken as —1 < a < 1 as usual. 

Using Boyer-Lindquist (Boyer & Lindquist 1967) co-ordinates, the corresponding metric element for the Kerr family of solutions in the 
spherical polar co-ordinate can be expressed as 



A /ir 



sin^ 6 



+ r''sme[l + — + — d0^ (12) 



where 6 is the polar angle, m = 1 + 77 cos^ 6 and A = — 2r -|- a^. 
The corresponding co-variant metric components are 



M 2 \ ^7-2 2 

fltt = - 1 • ffrr = — T-. gee = . 



gt4> = g<t>t 



2a sin^ 9 



2 . 2 . A , 2a^ sin'' 6 
= r sm I 1 + — H — 
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Corresponding contravariant elements can be expressed as 



tt 

a 



9t4, - 9ug,t,4, 



2r f_ a^^"! 



9 = 



1 _ A 

9rr Hr'^ 

1 1 

9ee = — = — 2 

900 (Ur^ 



2a 



9t^ - 9tt9<l»i> ^Ar 



914, ~ 9tt9<P4> A sin^ 9 

We, however, will be working on the flow solution on the equitorial plane (as defined by \d — 7r/2| ^ 1 in Novikov & Thome (1973), 
hence the corresponding line element for the equitorial slice can be expressed as 

2 2 2 2 2 

dSeg = {gtt)eqAt + {grr)eqAr + {gee)eqdO + 2{gt4,)eqdtd((> + {g4>4,)eqA<j) ; (13) 

where the subscript 'eg' implies that the corresponding values are evaluated at the equitorial plane. 
We thus obtain, 

(.gtt)eq = - (^1 - ^) , (14a) 

(fl'rr-)eq = (14b) 

{gee)eq = r•^ (14c) 

i9t4,)eq = {g4>t)eq = - — , (14d) 

r 

{g<t"t>)^i = ^; (i4e) 

where A = r* + r'^a'^ + 2ra^. The corresponding contravariant metric elements can thus be evaluated as, 

i.9t4>)iq - {gtt)eqig<l>cl>)eq Ar^ 



(g 



rr )eq 



{gee)eq = , \ = ^ (15c) 

{gt<l>)lq - {gtt)eq{g4,4,)eq Ar 

(g'^'^)eq = (94>4>)eq _ (l ~ r) (15e) 

{9t4>)iq- (9tt)^Ag4,4,)eq A 

Hereafter, however, we drop all the ' eg' subsripts for the sake of bravity. Any or g'^" will thus explicitly imply that the corresponding 
metric element has been evaluated on the equitorial plane. Using the cylindrical polar co-ordinate the corresponding metric on the equitorial 
plane can be expressed as, 

ds^ = giivAx^dx^ 

= -^dt' + ^dr' + ^ (d0 - Lodt) + dz^, (16) 

where z — r cos 6,lj — 2ar/A and Qzz = 1- 

For the metric element expressed using (r, 9, 4>)' g(r,6,(t>) = det((7^i/) = — r"', whereas for ds^ expressed using (r, 0, z),g(r,4,,z) = —t^- 
Main calculations presented in this work will be based on the line element as expressed in eq. (16). 
In this work, the polytropic equation of state of the following form 

p = Kp^' (17) 

is considered to describe the flow, where the polytropic index 7 (which is equal to the ratio of the specific heats, Cp and C„, respectively) of 
the accreting material is assumed to be constant throughout the fluid. A more realistic flow model would perhaps requies the implementation 
of a variable polytropic index having a functional dependence on the radial distance of the form 7 = 7(r) (Meliani et. al. 2004; Ryu & 
Chattopadhyay 2006; Migone & McKinney 2007; Mondal & Basu 2011; Mukhopadhyay & Dutta 2012). However, we have performed our 
calculations for a reasonable broad spectrum of 7 and thus believe that all astrophysically relevant polytropic indices are covered in our 
analysis. 
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The proportionality constant K in eq. (17) is related to the specific entropy of the accreting fluid provided no additional entropy 
generation takes place. Subjected to the condition that the Clapeyron equation of the form 

p=—pT, (18) 
fxrUp 

holds (in addition to the eq. (17)); kn, M and mp ~ tuh being the lotally measured flow temperature, the mean molecular weight, and the 
mass of the singly ionised hydrogen atom, respectively, the entropy per particle of an ensemble may be expressed as (Landau & Lifshitz 
1981) 

1 7 
■ log K H + constant; 



7 — 1 7 — 1 

where the constant depends on the chamical composition of the acccreting matter. K in eq. (17) can now be interpreted as a measure of the 
specific entropy of the accreting matter. 

The specific enthalpy h can thus be formulated as 

ft=^, (19) 
P 



where the energy density e includes the rest mass density and internal energy. Hence 
The adiabatic sound speed Cs is defined as 



e = p+^ (20) 
7-1 



= f I) (21) 

V / constant enthalpy 

The energy-momentum tensor of an ideal fluid can be expressed as 

r'"' = (e+p)wV+pfl^''. 

Vanishing of the four divergence of the energy momentum tensor provides the general relativistic version of the Euler equation i.e. 

T.^" = 0. (22) 

The continuity equation is obtained from 

ipv'^) = 0. (23) 

We have defined two killing vectors — 5f and (p'' = 5^ corresponding to stationarity and axisymmetry of the flow. 
We now contract eq. (22) with 0^ to obtain, 

<Pn [if^+P)v''v'"].^ + ii>^,p,^g^" = 0. 

But (jfp,!^ = due to axisymmetry, hence 
which further provides, 

5^4(e + p)vV]., = 0, (24) 
since 0'' = 5^. Since gij.\;i, = 0, eq. (24) can be written as 

[g^^{e + p)v>'v% = Q; 

from where we obtain 

[6,^,-].,. = 0. (25) 
From eq. (25) one thus infers (pfihv'^ = hv^ and hence hv^,, the angular momentum per baryon for the axisymmetry flow, is conserved. 

It can also be shown (in Fishbone & Moncrief (1976)) that the quantity v^v*, which is the angular momentum per unit inertial mass, is 
conserved for an iso-entropic flow and hence the world lines along which v^v* remains constant is a solution of the general relativistic Euler 
equation. 

In a similar way, we contract eq.(22) with 

ii.Tlt" = 0, 

to obtain that hvt is conserved. Hereafter, hvt will be interpreted as the specific energy of the flow and wiU be denoted be £ (scaled in unit 
of moc^ using the system of unit adopted in this work). 

For polytropic adiabatic accretion, £" is a first integral of motion along a streamline, and thus can be identified with the relativistic 
Bemouli's constant (Anderson 1989). 

The angular velocity can be defined in terms of the specific angular momentum A, where 

Vt 
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n = — = - + ^"'^ _r[2a + X{r - 2)] 



+ gt<i>>^ A. - 2a\r 

The normalisation condition v^Vy, = — 1 can be expressed as 

V Vt +V Vr + V^V^ = —1. 

In terms of the angular velocity O and the specific angular momentum A, one thus writes 

V*Vt + V^Vr + ^v^{ — Xvt) = —1 (26) 



Since = — ^^r^, one can rewrite eq. (26) as 

as (by converting v* — ^ vt) 

and finally, by using the definition of CI, since 



v*vt(l - AH - w^) = -1, 
(ff**-As*'*)(l-An-«V = -l 



1 (1 _ An - v^)v^ = 1, 



fft0 - 9ttg<i,4, 



and as it can be shown that u = — AO (where u is the advective velocity i.e. the radial velocity of the fluid measured in frame which 

co-rotates with the fluid such that v'^ = u/ -s/ QrrO- — u^), see Gammie & Popham (1998)), vt can be obtained as 



The mass conservation equation (the cointinuity equation as defined by eq. (23)) provides 

1 



iV^pv"),^ = 0, (28) 



where g = det(g,ii^). We multiply equation (28) with the co-variant volume element ^/^d*x to obtain 

dt and 84, are not relevant sue to the stationarity and axisymmetry. Also due to the assumption of non-existence of convection current along 
any non-equitorial direction, no non-zero terms involving (for spheriacal polar co-ordinate) or (for flow studied within the framework 
of cylindrical co-ordinate) may be obtained. One thus have 

dr{V^pv'')drded(l) = Q, (29) 

for accretion studied using the spherical polar co-ordinate and 

^ri^/^pv'')drdzd(|) = 0, (30) 

for accretion studied using the cylindrical co-ordinate. 

One can integrate eq. (29) for (/) = ^ 27r and 9 = —He — >■ Hg; ±Hg being the value of the polar co-ordinates above and below the 
equitorial plane, respectively, for a local flow half thickness H. 2H/r is cnstant for flow thickness in spherical polar co-ordinate for conical 
wedge shaped flow, (see e.g., Abramowicz & Zurek (1981); Gammie & Popham (1998); Lu, Yu & Young (1995); Lu et al. (1997a,b); Lu & 
Yuan (1998); Nag et. al. (2012) and references therein) and can integrate eq. (30) for z = —H^ — >■ (where ±i/z is the local half thickness 
of the flow) symmetrically over and below the equitorial plane for axisymmetric accretion studied using the cylindrical polar co-ordinate to 
obtain the conserved mass accretion rate M in the equitorial plane as 

p«M(r) = M; (31) 

where A{r) is the surface area through which the inward mass flux is estimated. For spherical symmetry, A{r) = AirHgr^ (for not so large 
value of 9) and for cylindrical symmetry, A{r) = AirHzr. 

In standard literature of accretion astrophysics, the local flow thickness for an inviscid axisymmetric flow can vary in three different 
ways, with different degrees of complexity, (see, e.g. Nag et. al. (2012)). A constant flow thickness is considered for simplest possible flow 
configuation where the disc height H is not a function of the radial distance (Abraham et al. 2006). In its next variant, the axisymmetric 
accretion can have a conical wedge shaped structure (Abramowicz & Zurek (1981), Lu (1985, 1986); Lu, Yu & Young (1995); Lu et al. 
(1997a,b); Lu & Yuan (1998); Lu & Gu (2004)) where H is directly proportional to the radial distance as H = A^r, and hence the 
geometric constant Ah is determined from the measure of the solid angle substended by the flow. For the hydrostatic equiUbrium in the 
vertical direction, (Matsumoto et al. (1984); Frank et al. (2002); Das & Czerny (2011) and references therein) expression for the local flow 
thickness can have a rather complex dependence on the radial distance and on the local speed of propagation of the acoustic perturbation 
embedded inside the accretion flow. In the present work, we consider the accretion flow to be in vertical equilibrium, and assiune that the 
flow has a radius-dependent local thickness with its central plane coinciding with the equitorial plane of the black hole. The equations (Barai, 
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Das, & Wiita 2004) of motion apply to the equitorial plane of the black hole, whereas the hydrodynamic flow varaibles are averaged over 
the half thickness of the disc H. We follow Abramowicz, Lanza & Pervival (1997) to derive the disc height for our flow configuration, and 
obtain the expression for the local flow thickness to be 



H{r) 



7 + 1 



(7- 



_{l-il + cl)}{X^vl-a%Vt-l)}_ 

Using the expression for vt as obtained from its expression in eq. (27) (using the expresions for g^^'s at equitorial plane as mentioned on 
Eqs. (14)) i.e., by writing vt as 



Vt 



Ar^A 



_(1-m2)(A2, 

H{r) can be obtained in terms of the advective velocity u. 



4aAryl + A2r2(4o2 -r^A)) 



5 CRITICAL POINT CONDITIONS 

The two first integrals of motion, the conserved specific energy £ of the flow and the mass accretion rate M can thus be expressed as 

(7-1) 



7 - (1 + cj) 



Ar^A 



A^ - 4AarA + \H'^{^d?- - r^A) 



(32) 



M = 4nA^ H{r)p- 



■VT^' '''' 

by integrating the stationary part of the energy momentum conservation equation and the continuity equation, respectively. The set of equa- 
tions (32 - 33) can not directly be solved simultaneously since it contains three unknown variables u, Cs and p, all of which are functions 

of the radial distance r. We would like to express p in terms of Cs and other related constant quantities. To accomplish such task, we 

1 1 

make a transformation H — M"f~-^ A - - i . Employing the corresponding equation for the sound speed as well as the equation of state, 

1 1 

S = M7 T--1 isT 7-1 can be expressed as 



47r A 2 Cs 



(7-1) 



7-(l+c2 



The entropy per particle a is related to K and 7 as Landau & Lifshitz (1994) 



H{r). 



(34) 



1 



■iogJs: + 



7 



+ constant 



(35) 



7-1 ° 7-1 

where the constant depends on the chemical composition of the accreting material. The above equation implies that K in some sense 
is a measure of the specific entropy of the accreting matter. Hence H may be interpreted as the measure of the total inward entropy flux 
associated with the accreting material and thus we label H to be the entropy accretion rate. It is worth mentioning here that the concept of the 
entropy accretion rate was first introduced in Abramowicz & Zurek (1981); Blaes (1987) for accretion under the influence of the Paczynski & 
Wiita Paczynski & Wiita (1980) pseudo-Schwarzschild Newtonian like black hole potential. H is thus conserved for the shock free polytropic 
accretion and increases discontinuously at the shock (if forms). Since £ and H remains conserved along a streamline, the spatial derivative 
(since we are dealing the stionary flow) of £ globally for shocked flow as well as the shock free flow and of H globally vanishes only for 
shock free flow. However, even if the shock forms, the spatial derivative of H valishes locally, and hence dE/dr = Q holds separately for the 
pre and the post shock flow, where the pre and the post shock flow implies the transonic accretion solution passing through the outermost and 
the innermost saddle type critical points, respectively. This point will further be clarified in the subsequent sections. 

The relationship between the space gradient of sound speed and that of the advective velocity can now be estabUshed by differentiating 
eq. (34) 



dcs 
dr 



(7- 



1 + 7 



4 



2m 



2 + mpa 



du 
dr 



(36) 



Differentiation of eq. (32) provides another relationship between dcs/dr and du/dr. We substitute dcs/dr as obtained in eq. (36) into that 
relationship and finally obtain 



2c^ 



(7+1) 



A 



1 2 
- + - 



4ip 



X 
2 



dr 



(1-W2) (7-|-l)(l-w2)u 



2ip 



(37) 
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where 



= X'vl - a' - 1) , Va = 1 - ^ ,<y = 2\^vt - a" , 



A- 

_ 1 dA _^ A dQ. _\ dr ' dr , ^^^^ 
^ A dr (1 - fiA) dr {g^^ + Xgtcf,) 

eq. (37) as well as eq. (36) can now be identified with a set of non-linear first order differential equations representing autonomous dynamical 

systems, and their integral solutions will provide phase trajectories on the radial Mach number M vs r plane. The regular critical point 
condition for these integral solution is obtained by simultaneously making the numerator and the denominator of eq. (37) vanish. The critical 
point condition may thus be expressed as 



'(i^ (7 + 1) 1/- 



2x6 



1/2 

, M|(r=ro) 



2r (r - 1) + 4A 



1/2 



(39) 



Following the aforementioned criteria, one obtains a 'smooth' or 'regular' critical point for which u,Cs and their space derivatievs are 
regular. For non dissipative inviscid flow, such a critical point may only be either of saddle type which allows transonic accretion solution 
to pass through it or of centre type through which no physical transonic solution can be constructed. Howevr, another kind of critical point 
can also be obtained for which D = does not ensure A' = 0, and one is left with an 'irregular' or 'singular' critical pointr.M and Cs are 
contineous but their derivatives diverge at such critical points. A singular critical point is obtained at the point of inflection of the homoclinic 
orbit for multi-transonic flow on the M — r phase plane. We will discuss this issue in greater detail in the subsequent sections while describing 
the procedure to obtain the multi-transonic shocked accretion phase topology. 

eq. (39) provides the critical point condition but not the location of the critical point(s). It is necessary to solve eq. (32) under the critical 
point condition for a set of initial boundary conditions as defined by [£, A, 7, a]. The value of Cs and u, as obtained from eq. (39), may be 
substituted at eq. (32) to obtain a complicated non-polynomial algebraic expression for r = rc,rc being the location of the critical point. A 
particular set of values of [£, A, 7, a] will then provide the numerical solution for such algebraic expression to obtain the exact value of Vc- 
It is thus important to find out the astrophycially relevant domain of numerical values corresponding to £, A, 7 and a. An appropriate choice 
for the four parameter initial boundary conditions are thus (Barai, Das, & Wiita 2004) [l ^ f ^ 2,0 < A$;4, 4/3^7^5/3, — l^a^l]. 
Although to be mentioned here that an upper limit for the Kerr parameter has been set to 0.998 in some works, see, e.g., Thome (1974). We, 
in our work, however, do not consider any such interaction of accreting material with the black hole itself which might allow the evolution 
of the mass and the spin of the hole as was considered in Thorne (1974) to arrive at the conclusion about such upper limit for the black hole 
spin. Hence we theoretically consider the domain for the Kerr parameter to be [— l^o^l]. 

The aforementioned four parameters may further be classified intro three different categories, according to the way they influence the 
characteristic properties of the stationary transonic aolutions. [f , A,7] characterizes the flow, and not the spacetime since the accretion is 
assumed to be non-self-gravitating. The Kerr parameter a exclusively determines the nature of the spacetime and hence can be thought of as 
some sort of 'inner boundary condition in qualitative sense since the effect of gravity is determined within the full general relativistic frame- 
work only up to several gravitational radii, beyond which it asymptotically follows the Newtonian description, [f , A] C [f , A, 7] determines 
the dynamical aspects of the flow, whereas 7 determines the thermodynamic properties. To follow a holistic approach, one needs to study the 
variation of the relevant features of the transonic accretion on all of these four parameters. 

Having equipped with the values of , A, 7, a], one can compute the location of the critical points by solving the algebraic equation 
as obtained by the substitution of eq. (39) in eq. (32). For convenience, the four dimensional hypersurface spanned by [f ,A,7, o] can 
be projected onto '^C-z different two dimensional or ''Ca different three dimensional parameter submanifolds to identify the regions of the 
parameter sub-space for which a multi-transonic (multi-critical solutions with stationary shock) accretion flow can be obtained. The accretion 
solution, as already mentioned, may be mono-critical or multi-critical with three critical points where two saddle type critical points are 
separated by a centre type critical point. The nature of a given critical point (whether it is of saddle type or a centre type) can be examined 
using certain eigenvalue equations and it can naturally be argued that any critical point associated with a stationary transonic solution will 
perforce has to be of saddle type (Chaudhury, Ray & Das 2006; Goswami, Khan, Ray & Das 2007). 

For multi-critical solution, the criteria for the accretion flow to have three critical points is associated with the value of the entropy 
accretion rate H (obtained from eq. (34)), will be denoted by H)i hereafter, evaluated for the solution passing through the innermost saddle 
type critical point is greater than the value of H evaluated for the solution passing through the outermost saddle type critical point, which 
will be denoted by Ho. The reverse situation, i.e.. So > S, provides the stationary configuration for which accretion solution connecting the 
infinity with the event horizon can have one critical point (the innermost one). For such a M — r phase portrait, the accretion solution through 
the outermost critical point is a part of the homoclinic orbit which can not be connected with the solution passing through the innermost 
critical point, and hence the stationary accretion for such configuration is essentially monocritical even if one obtains three formal solution 
of the critical point determining algebraic expression. 

If [£, A, 7, ^] represents the region of the four dimensional parameter space for which one obtains three critical points 

('mc' stands for 'multi-critical'), [£, A, 7, a]mca^l^' ^'T' '^Imc denotes the region embedded in [£, A, 7, o]^^ for which stationary accretion 
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solution can have three critical points - 'mca' being the acronym used for the phrase 'multi-critical accretion' . Hence it is the , A, 7, a\^^^ 
which we are interested in to identify the shocked multi-transonic flow. 

The space gradient for the advective flow velocity at the critical point is computed by solving the following quadratic equation 

fduV ^^Uu\ ^^^^^ 



\dr J , \ dr I , 

^ ' l(r=rc) ^ ^ l(r=rc) 

where the respective co-efficients, all evaluated at the critical point Tc, are obtained as 

(1 - u^) 7+1 7 + 1 



^ Cs (1 - i-i) „ v^Vta ^ 1 , 2A^ a 
o\ = —7- 2V' "2 = ——, 03 = 1 d4 = d2 



2 , uvtSs 



65 = 



u 1 — u _ 

3m^ -1 ^4 ^ (7 - 1 - c^) ^ _ (7 - 1 - eg) X , 8263x01 



u(l-v?) 1-82 ai(l-M2) ' ° 2ci 2(1-52)' 

r — 1 2 uwtx (4A^wt — a^) ip — vtcr"^ 

T\ = T 1 -— , T2 — 



A r Alp <jip 



_ C7T2X _ 1 2{r-iy 2 vtadx 

'^^ 4ip ' '^^ A A2 r2 4^/; dr ' 

2 r 2 r / 1 2\ , 2 1 XI 1 rfx 
^— - { (7 - 1 - c.) n + t;,c,r3} 3 J - 3 



2t;tM 

re 



(7 + 1)(1-m2) 



— (7 - 1 - Cs) + CsT3 



(41) 



Note, however, that all quantities defined in eq. (41) can finally be reduced to an algebraic expression in Fc with real coefficients that are 
functions of [£, A, 7, a]. Hence {du/dr)^^^^ is found to be an algebraic expression in rc with constant co efficients those are non linear 
functions of [£, A, 7, a]. Once Vc is known for a set of values of [£, A, 7, a], the critical slope, i.e., the space gradient for u at rc for the 
advective velocity can be computed as a pure number, which may either be a real (for transonic accretion solution to exist) or an imaginary 
(no transonic solution may be found) number. The critical advective velocity gradient for accretion solution may be computed as 

s) -l^^/f^i. <«) 

by taking the positive sign. The negative sign corresponds to the outflow/self-wind solution on which we would not like to concentrate in 
this work. The critical acoustic velocity gradient (dcs/dr)^^^^ can also be computed by substituting the value of {^)^_^ in eq. (36) and 
by evaluating other quantities in eq. (36) at rc- 

The values of the advective velocity u and the sound speed Cs evaluated at the critical point indicates that the Mach number at the 
critical point is not unity, and not even a constant as well. eq. (39) implies that the Mach number at the critical point is a function of the 
location of the critical point itself, and hence for any [E, A, 7, a]£[£, A, 7, a]^^^ one obtains three different Mach numbers corresponding to 
the three critical points for a multi-critical stationary solution. It is easy to show that Mc = (u/cs)^^ < 1 for all values of rc for a transonic 
flow, whether it is monocritical or multi-critical. Since a regular sonic point is identified with the radial distance (measured from the event 
horizon) where the transonic solution makes a continuous M < 1 — > M > 1 transition, the Mach number at the sonic point must be equal 
to unity. Hence the critical points and the sonic points are not topologically (as well as numerically) isomorphic. Such distinction between 
the critical and the sonic point is observed for polytropic accretion in the hydrostatic equilibrium in the vertical direction only and not for the 
polytropic flow with wedge shaped conical geometry or with constant thickness. This is a manifestation of the fact that the expression for the 
flow thickness (the disc height) only for accretion in vertical equilibrium is a function of the non constant sound speed and the expression 
for such disc height is obtained using a set of simplified assumptions, hence the dependence of the flow thickness on Cs is not exact. For 
isothermal accretion, critical points and the sonic points are identical even for flow in hydrostatic equilibrium in the vertical direction, let 
alone for the conical flow and for flow with constant thickness, since the isothermal sound speed is a constant and hence its space gradient 
does not contribute in the calculation of the space gradient of the dynamical velocity (from which the critical point conditions are derived). 
Since we concentrate on the polytropic accretion in the hydrostatic equilibrium in the vertical direction, we need to find out the sonic point 
by numerically integrating the flow equations. Two out of the three critical points in a multi-transonic accretion are of saddle type and the 
third one is the centre type. A physically acceptable transonic solution, however, can be constructed only through a saddle type critical point. 
No centre type critical point allows any transonic flow solution to pass through it. hence every saddle type critical point is accompanied by 
a sonic point r,, generally located at a radial distance (measured from the event horizon) smaller than the respective critical point rc- The 
criteria < rc is always satisfied since (u/cs)^^ < 1 whereas {u/cs)^^ = 1 and for a smooth transonic accretion, the Mach number 
anti-correlates with the radial distance measured from the event horizon. 
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Figure 1. Phase topology con'esponding to the multi-transonic shocked accretion and its associated wind branches obtained for a prograde flow characterized 
by [£ = 1.00001, A = 2.6, 7 = 2.6, a = 0.215]. The radial Mach number and the radial distance in logarithmic scale and in units of GMgu/c^ measured 
from the black hole event horizon has been plotted along the abscissa and the ordinate, respectively. ABBiCDE is the transonic accretion solution constructed 
thi'ough the saddle type outer critical point B (located at a distance 10815. 150 in units of GMgu / measured from the horizon) and the coiresponding outer 
sonic point Bi (located at a distance 8437.850 in units of GMgjj / measured from the horizon), FBG is the associated transonic self-wind solution passing 
through the outer critical point. The homoclinic orbit FUCiDiU consists of the transonic accretion and its associated self-wind solutions passing through the 
saddle type inner critical point I (located at a distance 6.999 in units of GMbh I (? measured from the horizon). The corresponding inner sonic point (not 
shown in the figure) is located at a distance 5.985 in units of GMbh / measured from the horizon. The centre type middle critical point shown by an 
asterisk (*) and marked by L located at a distance 17.213 in units of GMbh I <? measured from the horizon. Among two formally obtained shock transitions 
CCi and DDi, the shock formed between the outer sonic point and the middle critical point is found to be the stable one. ABBiCCiDiU thus represents the 
actual multi-transonic accretion flow connecting infinity with the black hole event horizon and contains the integral solutions constructed through the outer 
and the inner sonic points, respectively. 

6 CONSTRUCTION OF A TYPICAL MULTI-TRANSONIC PHASE TRAJECTORY 
6.1 The phase portrait 

We now demonstrate how to construct the integral multi-transonic solution with a stationary shock. In figure 1 we plot one such flow 
topology for \E = 1.00001, A = 2.6, 7 = 1.43, a = .215]. The radial Mach number has been plotted along the the Y axis and the equatorial 
radial distance r scaled scaled in units of GMbh /(? and measured from the event horizon, has been plotted along the X axis in logarithmic 
(logj^g) scale. Three regular critical points are obtained for such a multi-critical flow configuration - the outermost saddle type critical point 
r°"* located at r = 10815.150, the centre type middle critical point r™**^ located at r = 17.213, and the saddle type innermost critical point 
formed at r = 6.999. [7-°"*, r™"*, r*"] are computed by numerically solving the algebraic equation obtained by substituting the value 
of [u, Cs]^ (as defined in equation 39) in eq. (32). The value of u and Cs at any one of the three different critical points [r-°"\ r™"^, r*"] can 
now be calculated by substituting the respective values of the critical points back to the eq. (39). As already mentioned, Mach numbers can 
have different values computed at three different critical points. Hence r°"*, r™''*, and r™ are not co-linear on M — logj^Q(r) phase plane. 
The space gradient of the advective velocity (dit/dr)^. is computed using eq. (42) and the space gradient for the polytropic sound speed is 
calculated by substituting the value of {Au/dr)^ into eq. (36) in the critical limit. We now use the initial values [u, Cs, Au/dr, dCs/dr]j.out to 
integrate eq. (36) - (37) simultaneously to obtain the stationary transonic branch ABBiDE on the [M — log\o) phase plane passing through 
the outermost saddle type critical point r°"* as denoted by B on the phase plane. Since a critical point and a sonic point does not form at 
same radial distance, AB is not the subsonic branch of the flow, neither the segment BCDE represents the supersonic branch. The location 
of the sonic point is found by integrating du/dr and dcs/dr simultaneously to compute the radial equatorial distance for which Mach 
number becomes exactly equals to unity. For [£ = 1.00001, A = 2.6, 7 = 1.43, a — 0.215] used to obtain figure 1, we found the value of 
the outermost sonic point r°"* to be 8437.850 which has been identified in the figure 1 as Bi. Hence ABBi represents the subsonic flow 
and BiCDE represents the supersonic flow. We define Arcs ~ ('"c — '"s) to be a measure of the difference of the location of the critical 
and the sonic points. Hence the line segment measured along the X axis and corresponding to BBi represents the logarithmic value of 
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Arcs""* for the particular values of the [£,\, 7, a]used to obtain figure 1. One can study the dependence of Arcs""* on [£, A, 7, a]as well 
for the entire domain of the four parameters f , A, 7 and a to apprehend the effect of the black hole space time as well as the dynamical and 
the thermodynamic properties of the flow on the distinction between the sonic and the critical points. It is important to note that however 
small can Arcs be made, it never vanishes for any value of [£, A, 7, a]. This indicates that the non isomorphism of the critical and the sonic 
properties of the flow are not any artifact of the choice of the initial boundary conditions describing the stationary transonic accretion. 

If the value of the M is also provided along with [£, A, 7, a], one can calculate the values of all possible thermodynamic quantities 
corresponding to the flow, the pressure p, the density p and the ion temperature T of the accreting fluid, for example, at all radial distances 
stating from the infinity upto a very close proximity of r+ = 1 + Vl — o?- The transcritical solution passing through the outermost critical 
point r°"* seems to be doubly degenerate as is observed from the appearance of the phase topology FBG on the M — logj^(,(r) phase 
plane. FBG is obtained by integrating eq. (36 - 37) using [u, Cs, dtt/dr, dcs/dr] out but for the values of (dw/dr) out corresponding to the 

c c 

negative sign in eq. (42). Such twofold degeneracy is the consequence of the ±w degeneracy appearing in the expression for the energy first 
integral of motion as defined in eq. (32). Such degeneracy has, however, been apparently removed by orienting the phase portrait so that each 
phase topology represents either the accretion or the wind. The wind branch FBG, obtained by the advective velocity reversal symmetry, is 
a mathematical counter part of the accretion flow and is usually termed as the 'self -wind'. Had it been the situation that instead of starting 
from the infinity and heading toward the compact object, transcritical solution would generate from the close proximity of the accretor and 
would fly off from such object, FBG would denote the phase trajectory along which it would escape to infinity. The phrase 'wind solution' 
stems from the fact that the phase portrait corresponding to the solar wind solution due to Parker (1965) was topologically similar with the 
aforementioned mathematical counterpart of the accretion solution associated with the classical Bondi (1952) flow. 

Similar procedure may be used to obtain the transonic stationary accretion and the wind solutions passing through the innermost saddle 
type critical point rc" which is located at a radial distance r = 6.999 and is marked by I on the M — logio(r) phase plane. The corresponding 
value of the radial distance for the sonic point r'" comes out to r = 5.985. 

The transcritical accretion solution HIJ constructed through the inner critical point r™ folds back onto itself and joins with the corre- 
sponding transcritical self-wind branch HIK. The combined transcritical accretion-wind solution through the r'" thus forms a homoclinic 
orbit^ on the M — logio (r) phase plane. Such a homoclinic phase trajectory encompasses the centre type critical point r™"* flanked between 
rc"' and rl". The point of inflection H of the homoclinic orbit is actually a 'irregular' and 'singular' critical point. A tangent drawn through 
such a point of inflection comes out to be parallel to the Y axis. The Mach number M is defined at that point, whereas its space gradient 
dM/dr is not. The advective velocity u is continuous at that point but its space gradient diverges. At the point of inflection of the homoclinic 
orbit, the denominator of eq. (37) vanishes, allowing the corresponding numerator to assume a non-zero value. It is therefore understood 
that along with three regular critical points [r""*, r^*'', Tc"] , multi-critical stationary flow solution always possesses one more critical point 
which is of singular type. The only exception observed for a very special case where the multi-critical flow consists of two heteroclinic 
orbits'* since no homoclinic orbit for such configuration can further be realized. The transcritical heteroclinic orbits on M — logio{r) phase 
plane is characterized by the identical value of the entropy accretion rate H evaluated for the solution passing through the innermost saddle 
point r™ as well as for the solution constructed through the outermost saddle point r""*. 

A homoclinic orbit has its existence only in isolation and such trajectory does not qualify as a global transcritical solution since any 
realistic transcritical solution has to connect infinity with the event horizon to ensure the existence of the corresponding transonic flow. A 
local transcritical homoclinic integral flow solution can be made physically realizable by joining it with the transcritical non-homoclinic 
solution constructed through the outermost saddle type critical point r""' through a discontinuous shock transition since for non-dissipative 
inviscid flow two different transonic solution can not be smoothly connected to each other through any regular transition. In connection to 
the astrophysical flow. Such statement thus translates to the fact that no regular smooth stationary transonic solution can encounter more 
than one sonic point, and multi-transonic solution can only be realized when two different smooth transonic solutions can be connected 
through a stationary shock. The entropy accretion rate Hi„ for the accretion solution HIJ is greater than the entropy accretion rate Ho„t for 
the accretion solution ABiCDE subjected to the appropriate perturbative environment - a standing shock for which AH — — Sout^ 

amount of entropy accretion rate would be increased for the transonic solution passing through the outer sonic point Bi, would allow the 
aforementioned solution to make a discontinuous transition onto its subsonic homoclinic counterpart, i.e., the subsonic part of the transonic 
accretion solution HIJ. The combined multi-transonic shocked accretion solution would thus be consists of a segment (both subsonic and 
supersonic) of ABBiCDE and a segment (both subsonic and supersonic) of HIJ connected by a discontinuous shock, the location of which 
is to be determined by solving certain set of algebraic equations. 



A homoclinic orbit on a phase portrait is realized as an integral solution that re-connects a saddle type critical point to itself and embarrasses the corre- 
sponding centre type critical point. For a detail description of such phase trajectory from a dynamical systems point of view, see, e.g., Jordan & Smith (1999); 
Chicone (2006); Strogatz (2001). 

^ Heteroclinic orbits are the trajectories defined on a phase portrait which connects two different saddle type critical points. Integral solution configuration on 
phase portrait characterized by heterocUnic orbits are topologically unstable (Jordan & Smith 1999; Chicone 2006; Strogatz 2001). Subjected to the slightest 
possible perturbation, the heteroclnic loop opens up by forming a homoclinic orbit either through the inner saddle type point or through the outer saddle type 
point. 
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6.2 The relativistic shock of Rankine Hugoniot type 

In the present work, the first integrals of motion are the conserved specific energy and the mass accretion rate. For the non-dissipative inviscid 
accretion considered in our work, the shock produced is of energy preserving Rankine Hugoniot (Rankine 1870; Hugoniot 1887a,b; Landau 
& Lifshitz 1987; Salas 2006) type. The corresponding shock thickness has to be negligibly small compared to any characteristic length 
scale of the flow to allow no dissipation of energy as a consequence of the strong temperature gradient in between the inner and the outer 
boimdaries of the shock, where the terms 'inner' and the 'outer' are referred with respect to the proximity to the black hole event horizon. 
For a perfect fluid, we formulate the general relativistic Rankine Hugoniot condition to be 

[\puTu]\ = , 

[[TtMr7'']] = [[(p + £)«twr„]]=0, 

[[^^^'^n'^n]] = [[ip + e)w'r^ = o , (43) 

where = 1/ \/l — is the Lorentz factor and T^j^, is the corresponding energy momentum tensor. In the above equation, [[/]] denotes the 
discontinuity of any relevant physical quantity / across the surface of discontinuity, i.e., [[/]] = /2 — /i, where /i and /2 are the boundary 
values of the quantity / on the two sides of such surface. 

Simultaneous solution of eq. (43) yields the shock invariant quantity which changes continuously only across the shock surface. We 
obtain an analytical expression for such shock invariant quantity in terms of various local accretion variables and in terms of various initial 
bovmdary conditions describing the flow. During the numerical integration of the flow equations along the transonic solution ABBiCDE, 
we calculate the shock invariant. Simultaneously we calculate the same invariant while integrating the flow equations along the solution JIH 
starting from the inner sonic point upto the irregular sonic point on the homoclinic orbit (the point of inflection). We then determine the 
radial distance rgh where the numerical values of the shock invariant quantity is evaluated by integrating the two different flow segments as 
described above becomes identical. For every [£, X, 7, a]allowing the formation of a stationary shock, one in general obtains two different 
values of Vsh- Out of the two formal shock locations, the inner one (with reference to the proximity of r+) is always located in between the 
innermost and the middle sonic point, whereas the outer shock location is obtained in between the middle and the outermost sonic point. The 
shock strength M_ /M+ is different for these two shocks. Following the standard stability analysis procedure as provided in Yang & Kafatos 
(1995), one finds that the outer shock location is stable. Hereafter, we will refer to the stable outer shock location only whenever we use the 
phrase shock, and all shock related calculations will be performed with respect to that stable shock location only. 

ABBiCCiIJ on the M — logj^(,(r) phase plane as shown in figure 1. We need to calculate the values of various accretion variables along 
this segment of the flow topology with the information in mind that a sudden discontinuous transition for all such variables at the shock 
location has to be accounted for. We define the 'pre-shock variables' to be the value of any accretion variable at the shock location evaluated 
at the point C, and denote all such variables by a subscript Similarly, a 'post shock variable' is defined to be the value of the same variable 
(for same [£, A, 7, a]as well) evaluated at the point Ci and is denoted using a subscript '+'. The ratio of the pre (post) to the post (pre) shock 
variable for a fixed value of [f , A, 7, o]will provide the measure of the discontinuous change of such variables due to the presence of the 
shock. Had it been the situation that the shock would not form, and the flow would uninterruptedly follow ABB 1 CDF to approach the event 
horizon, the flow variables would also change continuously and no abrupt considerable alteration of the flow variables would be realized. 
Since a sudden change of the value of a flow variable is associated with the formation of a stationary shock in our model, a careful study 
of the radial profile of any accretion variable would provide a conclusive information about the appearance of a stationary shock. The radial 
variation of certain accretion variables (density, velocity, ion temperature etc.) are required to construct the observed spectra emergent from 
the accreting black hole, and the presence of shock can thus be inferred by investigating such spectral profile. 



6.3 Shock induced discontineous transition of flow variables 

In figure 2, we show the variation of the advective velocity u scaled in units of 10^*^ cm/sec, the bulk ion temperature T scaled in units of 10^° 
degree Kelvin, flow density p in gm/cc and fluid pressure p in dyne/cm^ for a black hole with mass Mbh = 3.6 x 10^ Mq and accretion 
rate M — XYZ. [£ = 1.00001, A — 2.6, 7 = 1.43, a = 0.215] has been used as the four initial boundary conditions to set up the flow. B 
and I are the outermost and the innermost saddle type critical points, respectively, and CCi indicates the shock transition. ABCE indicates 
the variation of the respective accretion variables along the transonic solution passing through the outer sonic point. If the shock would not 
form, the value of the respective variables would change continuously and monotonically, and the space evolution of such variables would be 
presented by the line segment ABCH. One could integrate the flow equations along the solutions passing through the outer sonic point upto 
the very close proximity of the event horizon and can obtain the value of the respective variable on ABCE ar an extreme close vicinity of the 
black hole event horizon for a shock free solution. If, however, the shock forms, there will be a sudden change of the value of the respective 
variable and its r variation profile can actually be demonstrated by the combined segment ABCCiU. Once again, one can integrate the set of 
differential and the algebraic equations governing the flow upto the very close proximity of r+ and the corresponding value of the respective 
variable at a radial distance nearly equal to r+ can be obtained for a shocked multi-transonic integral solution. 
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Figure 2. Variation of tlie advective flow velocity u (upper left panel), rest mass density p (upper right panel), flow ion temperature T (lower left panel) and 
pressure (lower right panel) corresponding to the flow topology presented in the figure 1 and for a black hole with mass 3.6 X 10*' Mq. The solid vertical line 
in each panel represents the shock transition and the labehng alphabets are in one to one correspondence with those used in figure 1, see text for further detail. 



7 DEPENDENCE OF THE SHOCK LOCATION AND THE SHOCK RELATED QUANTITIES ON [£ , A, 7, a] 

In this section we study the dependence of the shock location r^h and various pre and post shock values of the accretion variables on 
[£, A, 7, a]. To study such dependence on any particular parameter of the parameter set [£, A, 7, a], on specific flow energy £ for example, 
other three parameters A, 7 and a have to kept constant for the entire range of £ for which such dependence is studied. Whereas a wide 
range of choice for [f, A, 7, a] is available to produce a monotransonic accretion, only a limited non linear region of , A, 7, a] allows the 
existence of a shocked multi-transonic flow. A continuous range of all [£, A, 7, a] can not be used to construct such flow configuration since 
the Rankine Hugoniot condition is satisfied only for a small range of [£, A, 7, a]^^^^ G [iS, A, 7, a\^^^, where 'mcas' stands for 'multi-critical 
accretion with shock' and 'mca' indicates the 'multi-critical accretion' in general. [£, A, 7, a] C [£, A, 7, a]^^^ thus provides a true stationary 
multi-transonic accretion. We thus use various 'patches' of the region [£, A, 7, tsj^^^^j, to study the dependence of the shock related entities 
on [f , A,7,a]. 

One, however, needs to understand that such a choice of [f, A, 7, a] will indeed provide the generic profile for the aforementioned 
dependence. Consider a set of fixed values [Ai, 71, ai] to study the dependence of, say, the shock location Tsh on the available (for which the 
shock forms, subjected to the fixed set [Ai, 71, a\\) range of the specific energy starting from £min to £max- Such 'Vsh — £' profile can also 
be explored for any other fixed set, say [A2, 72, 0.2] for which the Rankine Hugoniot condition gets satisfied, only with the obvious difference 
that the numerical values corresponding to £min and £max associated with the flow described by [A2 , 72 , C12] will be different as compared to 
the values of £min and £max corresponding to the initial boundary conditions as defined by [Ai , 71 , ai]. Hence for any set of values [A, 7, a] 
for which the shock forms, the dependence of Tsh on £ can be studied. Similarly, the dependence of any shock related entity on any one of 
the initial boundary conditions [£, A, 7, a] can be studied for a fixed set of values of the rest of the initial boundary conditions for which a 
stationary multi-transonic accretion configuration can be realized. 

We observe that the shock location co-relates with the specific angular momentum A and anti-correlates with the specific energy £ and 
the polytropic index 7. Such trends are independent of the black hole spin parameter, and hence remain the same for the maximally rotating 
Kerr as well as for a non rotating Schwarzschild black hole. The aforementioned dependence does not explicitly provides any information 
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Figure 3. Vaiiation of the shock location r^i^ (plotted along the ordinate) with the Kerr parameter a (plotted along the abscissa) for the prograde accretion 
characterized by [S = 1.00001, 7 = 1.43] for three different values of the specific angular momentum A = 2.6 (uppermost panel), A = 2.17 (mid panel) 
and A = 2.01 (the lowermost panel). The shock location non-linearly con'elates with the black hole spin. 



about the dependence of the shock related quantities on the nature of the space time metric''. In this work, we are, however, mainly interested 
to study how the properties of the post shock flow at the close proximity of the event horizon are influenced by the spin parameter of the 
astrophysical black holes. Since such spin parameter determines the space time metric, our motivation is to study how the properties of 



^ We are dealing with a non self gravitating accretion hence no back reaction is considered and the metric is determined exclusively by the properties of the 
black hole itself. 
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Figure 4. Variation of of the sliock location r^j, (plotted along the ordinate) with the Kerr parameter a (plotted along the abscissa) for the retrograde accretion 
characterized by [8 = 1.00001, A = 3.3,7 = 1-4]. 

the transonic black hole accretion are determined by the nature of the black hole metric. Hence in the subsequent sections we study the 
dependence of the shock location as well as other shock related properties on the Kerr parameter in somewhat great detail. 

7.1 Dependence of the shock location on black hole spin 

As already argued, the entire range of the Kerr parameter, for the prograde as well as for the retrograde flow, can not be studied for any single 
fixed set , A, 7] since no such fixed set is available for which the Rankine Hugoniot conditions are satisfied for the entire range of the Kerr 
parameters [O^a^l], and [— IJSaJsO], for the prograde as well as for the retrograde flow, respectively. This can easily be shown by plotting 
the region [£, A, 7, a\^^^ embedded in the entire four dimensional hypersurface [£, A, 7, a] D , A, 7, a]^^^ to ensure that no single value 
of , A, 7] is available for which even a multi-critical solution, let alone a multi-transonic shocked solution, exists for the entire range of 
the Kerr parameter. We have chosen three different fixed sets [£ = 1.00001, A = 2.6, 7 = 1.43] , [£ = 1.00001, A = 2.17, 7 = 1.43] and 
[£ = 1.00001, A = 2.01, 7 = 1.43] to cover a significant range of the low to moderately high (0.2 <a< 0.5), high (0.85 <a< 0.925), 
and very high (0.9655 ^ a ^ 0.99) values of the Kerr parameters, respectively, to study the dependence of Vsh on the black hole spin for 
prograde flow. Such values of [£, A, 7] is chosen to maximize the available range of the Kerr parameter (for which the shock forms) for three 
different segments mentioned above, and to avail such range by minimally varying the initial configuration. For three different ranges, only 
the specific angular momentum A, that too by a considerably small amount, has been varied for each set of initial [£, A, 7], by keeping the 
subset [£,7] at its fixed value. 

As a representative value, we take [£ = 1.00001, A — 3.3, 7 — 1.4] to cover a reasonably large range of the black hole spin for which 
the shocked multi-transonoic accretion solution can be constructed for the retrograde flow. In figure 3, we plot the shock location rsh as a 
function of the Kerr parameter a for there prograde flow with fixed set of [£ = 1.00001, 7 = 1.43] and for three different values of A (as 
shown in the figure) for three different ranges of the Kerr parameters (as mentioned in the figure) for which the multi-transonic shocked 
solutions can be constructed. We observe that the shock location non-linearly co-relates with the black hole spin for the prograde flow. Hence 
we infer that for similar initial conditions describing the flow, the shock forms closest to the event horizon for a Schwarzschild type black hole, 
and furthest to the event horizon for an extremally rotating Kerr black hole. We find the same 'r^h — a' profile for the retrograde flow as well. 
In figure 4 the shock location is plotted against the black hole spin for the retrograde flow characterized by [£ = 1.00001, A — 3.3, 7 = 1.4] . 



7.2 Dependence of shock induced flow variables on black hole spin 

We would like to study how the characteristic dynamical and the thermodynamic features of the post shock flow are influenced by the black 
hole spin. One way of looking at this problem is to study the ratio of the pre (post) to the post (pre) shock values of various accretion 
variables, since such ratio is a marker of how the presence of a stationary shock introduces a sudden change in the value of the flow 
variables which in turn make a observable difference in the characteristic black hole spectra. For any flow variable /, /_ denotes its pre 
shock value evaluated at the shock location on the transonic solution passing through the outer sonic point and /+ denotes its post shock 
value evaluated at the shock location on the transonic solution constructed through the inner sonic point. For prograde flow characterized 
by [£ — 1.00001, A = 2.6, 7 = 1.43], in figure 5 we plot the ratio of the pre to the post shock Mach number (M^/M+), and post to the 
pre shock flow temperature (T+/T-), density (p+/p-) and pressure (p+/p-), respectively. (M-/M+) is termed as the shock strength as 
mentioned earlier and (p+/p-) is termed as the shock compression ratio. The shock strength anti-correlates with the shock location. The 
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Figure 5. Variation of the sliock strength (upper left panel), the shock compression ration (upper right panel), the ratio of the post to the pre shock temperature 
(lower left panel) and the ratio of the post to the pre shock pressure (lower right panel) with the Kerr parameter a (for each panel plotted along the abscissa) 
for the prograde flow characterized by [£ = 1.00001, A = 2.6, 7 = 1.43]. 

closer the shock forms to the event horizon, the higher the gravitational potential energy liberated resulting the formation of a stronger shock. 
A strong shock also compresses the flow by a considerable amount. As a result (since the shock as well as the flow under consideration 
is assumed to be energy preserving) the temperature and the pressure of the flow also increases. Hence [{p+/p-) , {T+/T-) , {p+/p-)] 
anti-correlates with the shock location and hence with the black hole spin angular momentum (the Kerr parameter a) 

Identical trend is observed for two other ranges of the black ole spin parameters explored (characterized by [£ — 1.00001, A = 2.17, 7 
= 1.43] and [£ — 1.00001, A = 2.01, 7 — 1.43]) in this work for the prograde flow, as well as for the retrograde flow characterized by 
[£ = 1.00001, A = 3.3, 7 = 1.4]. Such ' [{M-/M+) , {p-/p+) , {T+/T-) , (p+/p-)] - a' for the retrograde flow has been shown in fig- 
ure 6. Hence weakly rotating (low A), hot (high £), purely non relativistic (high 7) flow forms the strongest shock (the shock forms closest 
to the event horizon) for black hole with any value of the intrinsic spin angular momentum, i.e., the Kerr parameter a, and accretion onto a 
Schwarzschild black hole undergoes a strong shock transition compared to the case when matter with the same dynamical and thermodynamic 
properties accretes onto a Kerr hole. This summarizes how the shock formation phenomena and the properties of a multi-transonic flow get 
influenced by the space time metric (characterized by a), the dynamical (characterized by [£, A]) and the thermodynamic (characterized by 
7) properties of the accreting material. 



8 THE INFLUENCE OF BLACK HOLE SPIN ON QUASI-TERMINAL VALUES 

In this work, the numerical value of any accretion variable V evaluated at a very close proximity rg = r+ -\- 5 (r+ = 1 + yjl~~a?, and 5 
being a small number lying within the open interval < <5 << 1) of the event horizon is dubbed as the corresponding 'quasi terminal value' 
of V, and is distinguished by a subscript 5. The quasi-terminal value Vs is computed by integrating the flow equations (along the stationary 
transonic branch) from the critical point Tc down to ra. We take 5 = Q.QQIGMbh /(? and perform our calculation for Vs for a 3 x \(f Mq 
black hole accreting at a rate of 4.29 x lO-e^A^eYr"^. Such values of Mbh and M corresponds to our Galactic centre black hole and 
its environment where low angular momentum inviscid advective accretion model is considered as an appropriate approximation (see, e.g., 
Moscibrodzka, Das & Czerny (2006) and references therein). Also to be noted that Mbh and M used in this work are two representative 
values, and any other value for the black hole mass as well as for the accretion rate can be considered for our calculation of Vs. 

For any generic flow variable V we calculate the corresponding Vs along two branches, either along the solution ABBiCCiII for a 
multi-transonic shocked accretion flow, or for a hypothetical shock free solution ABBiCDE passing through the outer sonic point only. For 
the same set of initial boundary condition [£, A, 7, a], the dependence y^^°'^'' and Vj"" ^^""^^ on the black hole spin angular momentum as 
well as on [£,\, 7] can be studied to estimate the impact of the shock formation phenomena in determining the properties of the matter 
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Figure 6. Variation of the sliock strength (upper left panel), the shock compression ration (upper right panel), the ratio of the post to the pre shock temperature 
(lower left panel) and the ratio of the post to the pre shock pressure (lower right panel) with the Kerr parameter a (for each panel plotted along the abscissa) 
for the retrograde flow characterized by [E = 1.00001, A = 3.3, 7 = 1.4]. 



extremely close to the black hole. This in turn will allow us to infer the influence of the shock formation on the observable spectra generated 
by the photon flux emanating out from the region inside the ISCO. 

For CO rotating flow, in figure 7 we plot the variation of the quasi-terminal values of the Mach number (Ms, the top left panel), 
flow density {ps in CGS units, top right comer), bulk ion temperature {Ts, in units of degree Kelvin, bottom left panel) and pressure 
{ps in CGS unit, bottom right panel) respectively, for both shocked solution (solid red line) and for the shock free solution (dashed blue 
line) on the black hole spin for multi-transonic flow characterized by [£ = 1.00001, A = 2.6,7 = 1.43]. Similar relational dependence 
are presented in figure 8 and figure 9 for other ranges of the Kerr parameter for which multi-transonic shocked flow can be described 
by [£ = 1.00001, A = 2.17, 7 = 1.43] and [£ = 1.00001, A = 2.01, 7 = 1.43], respectively. For the same set of [£, A, 7, a], we find that 
j^^shock ^ j^no shock pshock ^ ^no shock ^ whcreas Ti''""'' > rr and pf°'='' > pT 

From figure 1, one observes that for any r < r"\ the value of Mach number evaluated along the transonic solution passing through 
the outer sonic point is always greater than that evaluated on the transonic solution passing through the inner sonic point. At the shock, 
the Mach number decreases discontinuously, and although the value of the Mach number shoots up at a very high rate (space gradient of 
the Mach number, i.e., dhl/dr, becomes large), the Mach number for the supersonic flow in the post shock region can never exceed the 
value of the Mach number associated with the supersonic segment of the shock free solution at any r since in that case the post shock 
supersonic flow would have to intersect the shock free supersonic solution on 'M — logu){r)' plane. Such a crossover is not allowed since 
no two phase topology can intersect on a phase plane. Hence what actually one observes is [dM / drY'^°'^^ > (dM / dr)"° "^"'^'^ but the 
trend Mg'^""'' < Mg° "'^""'^ is maintained. At the extreme close proximity of the event horizon (rs << 0.001) the post shock supersonic 
branch asymptotically approaches the shock free supersonic branch, hence j\f|'*°='= — >Mg° for such an extremely small value of rg. 
Nevertheless, the criteria Mt""'' - '''""^VO remains valid for for all values of rg, however small rg can be made. 

Similar situation is observed for the variation of pg with the spin as well. Although the density increases at the shock, close to r+ the 
flow density corresponding to the shocked flow makes a crossover with the density profile corresponding to the shock free transonic flow, 
and starts decreasing gradually as has been observed in figure 2. For the flow temperature and pressure, no such crossover takes place and 
hence the trend Tt°-''' > and pf""'' > are uniformly maintainedi, see, e.g., figure 2. Similar ' [Mg, pg, Tg,pg] - a' 

profile is observed for the retrograde flow as well, see, e.g., figure 10 for such variations for the counter rotating accretion. 




Figure 7. Variation of tlie quasi-terminal Macli number (upper left panel), quasi-terminal density (upper right panel), quasi-terminal temperature (lower left 
panel) and the quasi-terminal pressure (lower left panel) with the KeiT parameter a (plotted along the abscissa) for shocked (dashed red line) and for the 
hypothetical shock free (solid black line) prograde flow characterized by [£ = 1.00001, A = 2.6, 7 = 1.43]. 




Figure 8. Variation of the quasi-terminal Mach number (upper left panel), quasi-terminal density (upper right panel), quasi-terminal temperature (lower left 
panel) and the quasi-terminal pressure (lower left panel) with the Kerr parameter a (plotted along the abscissa) for shocked (dashed red line) and for the 
hypothetical shock free (solid black line) prograde flow characterized by [£ = 1.00001, A = 2.17, 7 = 1.43]. 
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Figure 9. Variation of tlie quasi-terminal Macli number (upper left panel), quasi-terminal density (upper right panel), quasi-terminal temperature (lower left 
panel) and the quasi-terminal pressure (lower left panel) with the KeiT parameter a (plotted along the abscissa) for shocked (dashed red line) and for the 
hypothetical shock free (solid black line) prograde flow characterized by [S = 1.00001, A = 2.01, 7 = 1.43]. 




Figure 10. Variation of the quasi-terminal Mach number (upper left panel), quasi-terminal density (upper right panel), quasi-terminal temperature (lower left 
panel) and the quasi-terminal pressure (lower left panel) with the Kerr parameter a (plotted along the abscissa) for shocked (dashed red line) and for the 
hypothetical shock free (sohd black line) retrograde flow characterized by [S = 1.00001, A = 3.3, 7 = 1.4]. 
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Figure 11. Variation of the quasi-terminal Mach number (upper left panel), quasi-terminal density (upper right panel), quasi-terminal temperature (lower left 
panel) and the quasi-terminal pressure (lower left panel) with the specific flow energy E (plotted along the abscissa) for shocked (dashed red line) and for the 
hypothetical shock free (solid black line) prograde flow characterized by [A = 2.17, 7 = 1.43, a = 0.881049812] . 



8.1 Dependence of quasi-terminal values on , A, 7] 

We also study the dependence of the quasi-teiminal values on the specific energy E (figure 11), specific angular momentum A (figure 12), 
and the flow adiabatic index 7 (figure 13). For all such cases the following trend is found 

M shock , tl rno shock shock , no shock mshock ^ rr^no shock shock ^ no shock / a a\ 

6 < Ms ,ps < Ps ,Ts >Ts ,ps > P6 (44) 

as has been observed for the black hole spin dependence of [Ms, ps,Ts,ps]- 



9 SPIN DEPENDENCE OF QUASI-TERMINAL VALUES FOR THE STATIONARY MONO-TRANSONIC ACCRETION 

It is instructive to investigate whether the characteristic feature of the black hole spin dependence of the quasi-terminal values remain 
invariant for a direct spin flip of the astrophysical black hole. In other words, we would like to understand whether the ' [Ms, ps, Ts,ps\ — a' 
profile gets significantly altered when the initial boundary conditions are switched from [f, A, 7, a] to [£, A, 7, — a]. One understands that the 
aforementioned exercises can not be performed for the multi-transonic accretion since the Rankine Hugoniot conditions can not be satisfied 
for a certain [£, A, 7, a]as well as for the same (magnitude wise) values of the [£, A, 7] and a but with the directionally opposite value of a, 
i.e., for [£, A, 7, —a]. This is a consequence of the fact that [£, A, 7, does not allow any parameter degeneracy. 

Construction of the stationary mono-transonic solutions, however, are not constrained by such limitations and one can obtain such 
solutions when [£, A, 7, ajgets directly swapped to [£, A, 7, —a]. We choose a suitable set of values of [£ = 1.2, A = 2.0, 7 = 1.6] for which 
a stationary transonic accretion solution can be constructed for the entire range of the black hole spin, i.e., for — l^a^l. A considerably 
slowly rotating substantially hot and almost purely non relativistic flow allows us to study the black hole spin dependence of the quasi-terminal 
values for the entire range of the Kerr parameters describing both the prograde and the retrograde accretion. For such flow configurations, we 
obtain that the critical as well as the sonic points are always of the innermost type saddle one. The location of the critical and the sonic points 
are found to be at close proximity of r+ for the prograde accretion onto a maximally rotating hole, and are formed at the maximally allowed 
distance of approximately (4 — 4.5) GMbh / (? unit for the retrograde accretion onto maximally rotating hole with the negative value of the 
spin parameter. The value of Ar^ = {tc — Ts) thus comes out to be minimum for a — > 1 and maximum for a — > —1. For mono-transonic 
accretion, in figure 14 we show the dependence of [Ms, pdslta, Ts,ps] on the black hole spin parameter for the entire range of — l^a^l. 
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Figure 12. Vaiiation of the quasi-terminal Mach number (upper left panel), quasi-terminal density (upper right panel), quasi-terminal temperature (lower left 
panel) and the quasi-terminal pressure (lower left panel) with the specific flow angular momentum A (plotted along the abscissa) for shocked (dashed red line) 
and for the hypothetical shock free (solid black line) prograde flow characterized by [€ = 1.000004, 7 = 1.43, a = 0.881049812] . 
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Figure 13. Variation of the quasi-terminal Mach number (upper left panel), quasi-terminal density (upper right panel), quasi-terminal temperature (lower left 
panel) and the quasi-terminal pressure (lower left panel) with the adiabatic index 7 (plotted along the abscissa) for shocked (dashed red line) and for the 
hypothetical shock free (soUd black line) prograde flow characterized by [£ = E = 1.000004, A = 2.28, a = 0.881049812]. 
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Figure 14. Variation of the quasi-terminal Mach number (upper left panel), quasi-terminal density (upper right panel), quasi-terminal temperature (lower 
left panel) and the quasi-terminal pressure (lower left panel) with Kerr parameter a (plotted along the abscissa) for mono-transonic flow characterized by 
[£ = 1.2, A = 2.0,7 = 1-6]. The entire range of the Kerr parameter for both the prograde as well as the retrograde flow has been taken care of. The prograde 
and the retrograde branches are not symmetric. 



10 LINEAR STABILITY ANALYSIS OF THE STATIONARY SOLUTION 

In this work, our entire analysis of the multitransonic flow at the viscinity of the hole is based on the phase space behaviour of the stationary 
solutions. It is, however, important to ensure that such stationary configurations are stable as well, at least upto the limit of astrophysically 
relvant time scales. This can be accomplished by studying the time evolution of a linear acoustic like perturbation applied around a stationary 
configuration. We will first demonstrate that the axisymmctric accretion can be considered as potential flow and will thus identify the 
corresponding velocity potential. Next we perturb such velocity potential and will prove that such perturbation will not diverge to destabilize 
the original stationary solution of our interest. In the folowing part we present the analysis for irrotational, entropy conserving flow solutions 
in general, which necessarily includes the stationary ones. 
From Thermodynamics, 

dfe = TdQ+^^ (45) 

For polytropic flow along a specified streamline, we have from eq. (22) , 

hv^ipv"):,, + pv^ihv").,, + p^^g"'' = 0. (46) 

Since due to speecific entropy conservation along a streamline, dh = from eq. (23) one obtains, 

v''{hv^);,. + d^h = 0. (47) 

We define Qai3 = PUP^'jJij.v to be the vorticity of the flow, where uj^,, = (hv/j,);^ — {hv^)-,^, and = 5^ + VaV^ is the projection tensor. 
Since Q,ap = for irrotational flow, we obtain 

0Ja,3 + V[jv'^LJav + VaV^UJf_,fj + VaV'^ V 0v'^ UJ = 0. (48) 

The 2nd, 3rd and the 4th terms in the above expression vanish owing to the relation w^w^ = — 1 and by virtue of eq. (47). We thus obtain 
ujai3 = which implies 

{hvc);p - {hvp);c = 0. (49) 

eq. (49) indicates that there exists a 4-scaIar * such that, 

hVc = -^a*. (50) 
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^ is thus the velocity potential of the flow we consider in this work. 

Solutions (subjected to different initial boundary conditions) describe the vorticity free polytropic flow. Evidently the flow is entirely 
determined by the solutions corresponding to the velocity potential ^ . To analyse the stability of such solutions we may introduce small 
perturbations on the solutions which satisfy the eq. (50) and investigate whether such perturbations may eventually grow to mask the original 
solutions. If such perturbation does not diverge, the stationry solutions are proved to be stable within our framework. 

Applying small perturbation on the backgroimd values of the flow variables as, /i — >■ /i + 5h, p ^ p + Sp, V/j, ^ V/j, + Sv^i, we obtain, 

to^ + 5hv^ + hSv^ = 5^* + a^(5*. 

Perturbation of the eq. (23) provides, 

{5p v^),^^ + ipSvn,^ = 0. 
Using the relation, v'^Svi_i = 0, and denoting by /, we obtain, 

1 



da 



1 - ^ )■ v'^v" 



: 0. 



(51) 



For axisynmietric flow configuration as considered in our work, 

v'^ = (v^,v^ ,0, v'^^ 

and from radial propagation of perturbation we obtain, 

dj= (dtf,drf,0,0) 

Substituting these conditions into the eq. (51), it reduces to, 

d^.ir^d.f) = 0; (52) 

where p and u run for and 1 and 

eq. (52) is the equation which determines the time evolution of the first order linearly perturbed velocity potential (see Moncrief (1980); BiUc 
(1999)). 

Substitution of the trial acoustic wave solution of the form / = f^' cxp(— twt) into the equation (52) yields. 

The space dependent part is expressed in terms of the trial power series of the following form, 

(r) = exp 2^ 



and it is examined whether the solution is bounded within the finite limit at the outer boundary as r — ^ oo. Collecting the coefficients of the 
same power of w (w >> 1) together we obtain. 



for w^containing terms, - /" - 2it^^ + f ( ' = 0, 

dr \ ar J 

, 1 , . . , „.,t.dfco .df^ ^ drdk-i 

tor uj containing terms, — 2tf — i — 1 h / 

dr dr dr dr 



^ dfc_i dfcp _^ d fc-i 



f ^ • • ^ o • ftr dfci , df'' dko r 

for w containing terms, — 2if — 1 — h / 

dr dr dr 



^ dfc-i dfci 
dr dr 



dr dr ' dr^ 
+ 



dfcp 
dr 



^ d^fco 



dr^ 



0, 
= 0. 



The leading order coefficients turn out to be (from eq. (53)), 

fc-i = 

and substituting back into the eq. (54) we obtain. 



I dr 

J /'•'■ 



(53) 
(54) 

(55) 
(56) 



dko 
dr 



dr dr 



which, upon further simplification becomes 



dfcp 
dr 



'- 2i Pr _ 2i frr f^^ ±^ ^^p- - frr 



frrft 



(57) 
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Hence we find, 



From eq. (55) we obtain 



ko 



dfci 
dr 



{Mil 

[ dr 



ifco I rrr f dfcp \ , rrr d^fcp 
dr ' J \ dr I 'J dr'^ J 



2[i/*'--/'-'^%i] 

From eq. (52) using expressions of contravariant metric elements as defined on the equitorial plane in Eqs. (15), one obtains, 

r = 
f = 



— p A A, 1 , 



(58) 
(59) 

(60a) 
(60b) 
(60c) 



It is easy to see that in the asymptotic limit r — >■ oo, A ~ and A ~ r''. At that limit p tends to a constant ambient value, denoted by 



3 . The Lorentz factor 71, tends to unity for 



Poo and subsequently h tends to its ambient value h^o, as also tends to some ambient value Cs 
accretion at the outer boundary condition. 

To find out the the asymptotic behaviour of u, we make use of the eq. (31). It turns out that 

u 

VI ~ pA' 

where A is 4:nHzr for the flow considered within the framework of the cylindrical symmetry. It is to be noted that is constant for the 
accretion with constant flow thickness. For axisymmetric accretion in vertical equilibrium, 

(7-l)c? 



7+1 



_ [7 - (1 + ci)] [X^gyh^ - a?{£/h - 1)]_ 

On the other hand for conical model described within the framework of spherical synmietry for —Hg ^ 9 ^ Hg, the area A is AirHgr^, 
where Hg is constant for a disc. Hence for all the flow configurations at asymptotic limit r — ^ 00, ^ ~ r" where a ^ 1. Thus, 



at this asymptotic limit. 

Now it becomes apparent from equations (60) that. 



f 

r 
r 



'-9 



(61a) 
(61b) 

(61c) 



Hence substituting these expressions (61) into eqs. (56), (58) and (59) one readily obtains that fe-i ~ r, fco ~ Inr and fci ~ 1/r at the 
asymptotic limit. Hence for the first three terms in the trial power series solution for the space dependent part of the perturbation it is evident 
thata;|if-i| >> |fc()| >> |fci |/cj for the high frequency regime at the asymptotic solution. This indicates that w"|A;„| >> 
or in other word the power series converges even at the outer bounadary of the flow ruling out any divergence of any possible perturbing 
component of 3* from the solution of eq. (50) for a particular set of boundary conditions. 

In this work, we study how the stationary accretion solutions in the Kerr metric behaves close to the black hole hole horizon. The 
dependence of such behaviour on the black hole spin is also studied. By employing suitable stability analysis above we ensure that such 
stationary solutions (which constitute a subcategory of the potential flow, in general) are stable, and hence any spectral profile which might 
be constructed out of those solutions are reliable, at least upto an astrophysicaUy relevant time scale. 



11 DISCUSSION 

Majority of the works on the role of the black hole spin in influencing the accretion dynamics are focused on high angular momentum disk- 
like flows, with the central role payed by the Innermost Stable Circular Orbit (ISCO). However, SgrA* and M87 are the two most preferred 
candidates for direct imaging of the flow close to the black hole horizon since the angular size of the black hole horizon is by far the largest in 
these two sources due to the interplay between the black hole mass and the distance to the source is measured by the observer. The exact value 
of the angular momentum of the inflowing material for the aforementioned two sources is difficult to estimate, and the current evaluations 
indicate values ranging from moderate (Cuadra et. al. 2008) to considerably low (Moscibrodzka, Das & Czemy 2006). 
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The present constraints inferred on tlie unresolved components in these two sources in the mm bands are already very impressive 
(Doeleman et al. 2008, 2012). In the future VLBI Event Horizon Telescope will bring us still much closer to the central compact object for 
these sources (Doeleman 2010). In order to understand the salient features of these images as well as the properties of the corresponding 
broad band radiation spectra we need to predict the emissivity distribution and construct the expected black hole silhouettes for various 
models using the ray tracing techniques as described in (Lvmunet 1979; Fukue & Yokoyama 1988; Karas et. al. 1992; Falcke, Melia & Agol 
2000; Takahashi 2004; Huang, Cai, Shen & Yuan 2007). 

Sizable amount of work in this direction has been performed for high angular momentum flows and for ion tori (e.g. Vincent et al. 
(2012)). The results given in our present paper paper form a starting point for complementary study for the case of low angular momentum 
accretion. In our next work (Huang and Das, in preparation), we plan to perform the black hole shadow imaging corresponding to the low 
angular momentum axisymmetric accretion as considered in this work. 

In our present work we have studied the radial ion temperature profile as a function of the Kerr parameter, as well as the dependence of 
the corresponding Tf, profile on the black hole spin. Our ongoing calculation concentrates on the calculation of the electron temperature from 
such ion temperature. Knowledge of such electron temperature, along with the density and the velocity profile as calculated in our present 
work, will then provide us the complete knowledge of the emitted polarized radiation in the millimeter and sub- millimeter band. We also 
plan to study the influence of shock formation on the polarized emission of SgrA* using our flow model. 

Findings as illustrated in Fig. 14 indicate a consistent asymmetry between the prograde and the retrograde flow as far as the quasi- 
terminal values are concerned. The constructed shadow image is also expected to manifest the asymmetry. We thus expect to propose a novel 
method to differentiate the co-rotating and counter-rotating flow once we construct the corresponding spectra and the associated shadow 
images out of the accretion variables as calculated using our flow model. The observations of Sgr A* and M87 by Doeleman et al. (2008, 
2012) seem to rule out counter-rotating flow since the images are in both cases much smaller than ISCO. However, the argument applies only 
withing the frame of disk-like accretion. Low angular momentum flow is much less influenced by the position of the ISCO and the emissivity 
is more concentrated towards the center as perceived in a spherically symmetric flow. Therefore, the construction of the predicted images in 
the case of low angular momentum flows requires urgent attention as we believe. 



ACKNOWLEDGMENTS 

SH and IM would Uke to acknowledge the kind hospitality provided by HRI, Allahabad, India, under a visiting student research programme. 
The visits of PB, SN and TN at HRI was partially supported by astrophysics project under the Xlth plan at HRI. The work of TKD has been 
partially supported by a research grant provided by S. N. Hose National Centre for Basic Sciences, Kolkata, India, under a guest scientist 
(long term sabbatical visiting professor) research programme, as well as is partially funded by the astrophysics project under the XI th plan 
at HRI. PB acknowledges support from the ERG Starting Grant "cosmoIGM". 



REFERENCES 

Anderson, M. 1989, MNRAS, 239, 19 

Abraham, H., Bilic, N., Das, T. K., 2006, Classical and Quantum Gravity, 23, 2371 
Abramowicz, M. A. & Zurek, W.H., 1981, ApJ, 246, 314 

Abramowicz, M. A., Czerny, B., Lasota, J. P., Szuszkiewicz, E., 1988, ApJ, 332, 646 

Abramowicz, M. A., & Chakrabarti, S. K., 1990, ApJ, 350, 281 

Abramowicz, M. A., Chen, X. M., Granath, M., & Lasota, J. P., 1996, ApJ, 471, 762 

Abramowicz, M. A., Lanza, A., & Percival, M. J. 1997, ApJ, 479, 179 

Abramowicz, M. A., Kato, S. 1989, ApJ, 336, 304 

Abramowicz, M. A., Zurek, W. H. 1981, ApJ, 246, 314 

Afshordi, N., Paczyhski, B., 2003, ApJ, 592, 354 

Artemova, I. V., Bjomsson, G., Novikov, I. D., 1996, ApJ, 461, 565 

Armitage, R J., Reynolds, C. S., & Chiang, J., 2001, APJ, 648, 868. 

Barai, P, Das, T. K., Wiita, P J., 2004, ApJ, 613, L49 

Bardeen, J., M., Press, W., H., & Teukolsky, S., A., 1972, ApJ, 178, 347 

Becker, R A., & Le, T., 2003, ApJ, 588, 408 

Beloborodov, A. M., & Illarionov, A. F 1991, MNRAS, 323, 167 

Bhattacharjee, J. K., Ray, A. K., 2007, ApJ, 668, 409 

Bhattacharya, S., Bhattacharjee, J. K., 2005, Proceedings of the Indian National Science Academy, 71A, 1 

Bihc, N., 1999, Class. Quantum Grav., 16, 3953 

Bisikalo, A. A., Boyarchuk, V. M., Chechetkin, V. M., Kuznetsov, O. A., & Molteni, D. 1998, MNRAS, 300, 39 
Blaes, O., 1987, MNRAS, 227, 975 

Bohr, T., Dimon, P., Putkaradze, V., 1993, J. Fluid Mech., 254, 635 



Accretion close to event horizon 29 



Bondi, H. 1952, MNRAS, 112, 195 

Boyer, R. H., & Lindquist, R. W. 1967, J. Math. Phys. 8, 265 

Bray, A. J., & Moore, M. A. 1987, Phys. rev. Lett. 58, 57 

BuUga, S. D., Globina, V. I., Gnedin, Y. N., Natsvhshvih, T. M., Pitrovich, M. Y., & Shakht, N. A., 2011, Astrophysics, Volume 54, Issue 
4, pp.548 

Caditz, D. M., & Tsuruta, S. 1998, ApJ, 501, 242 
Chakrabarti, S. K., 1989, ApJ, 347, 365 
Chakrabarti, S. K., 1990, ApJ, 350, 275 
Chakrabarti, S. K., 1996, MNRAS, 283, 325 
Chakrabarti, S. K., 1996, ApJ, 471, 237 

Chakrabarti, S. K., 1990, Theory of Transonic Astrophysical Flows, World Scientific, Singapore 

Chakrabarti, S. K., Titarchuk, L. G., 1995, ApJ, 455, 623 

Chattopadhyay, I., & Das, S., 2007, New Astronomy, Volume 12, Issue 6, p. 454 

Chaudhury, S., Ray, A. K., Das, T. K., 2006, MNRAS, 373, 146 

Chen, X., Abramowicz, M. A., Lasota, J. P, 1993, ApJ, 476, 61 

Chen, X., Taam, R., 1993, ApJ, 412, 254 

Chicone, C, 2006, 'Ordinary Differential Equations with Applications', Springer; 2nd edition. 

Cuadra, J., Nayakshin, S., Martins, R, 2008, MNRAS, 383, 45 

Czerny, B., Lachowicz, P, Dovciak, M., Karas, V, Pechacek, T, & Das, T. K., 2010, A & A Volume 524, id.A26 
Daly, R. A., 2011, MNRAS, 414, 1253 

Das, S., 2007, MNRAS, 376, 1659 

Das, S., & Chattopadhyay, I., 2008, New Astronomy, Volume 13, Issue 8, p. 549 

Das, T. K., & Chakrabarti, S. K. 1999, Class. Quantum Grav. 16, 3879 

Das, T. K., 2002, ApJ, 577, 880 

Das, T. K., & Sarkar, A., 2001, A & A, 374, 1150 

Das, T. K., Pendharkar, J. K., Mitra, S., 2003, ApJ, 592, 1078 

Das, T. K., Rao, A. R., & Vadawale, S. R. 2003, MNRAS, 343, 443 

Das, T. K., & Czerny, B., 2011, MNRAS, 414, 627 

Das, T. K., & Czerny, B., 2012a, MNRAS, 421, L24 

Das, T. K., & Czerny B., 2012b, New Astronomy 17, 254 

Das, T. K., Huang, L. M., & Inoue, M., 'Black hole shadow imaging for shocked accretion flow in Kerr metric'. (In Preparation). 
Dauser, T., Wilms, J., Reynolds, C. S., & Brenneman, L. W, 2010, MNRAS, 409, 1534 

Doeleman, S.S., et al. 2008, Nature, 455, pp. 78 

Doeleman, S., 2010, "Proceedings of the 10th European VLBI Network Symposium and EVN Users Meeting: VLBI and the new generation 
of radio arrays. September 20-24, 2010. Manchester, UK. PubUshed online at http://pos.sissa.it/cgi-bin/reader/conf cgi?confid=125, id.53" 
Doeleman, S.S., et al., 2012, Science, 338, pp. 355 
Falcke, H., Meha, R, & Agol, E., 2000, ApJ 528, L13 
Ferrari, A., Trussoni, E., Rosner, R., and Tsinganos, K. 1985, ApJ, 294, 397 
Fishbone, L.G. & Moncrief V, 1976, ApJ, Volume 207, pp. 962 

Frank, J., King, A., Raine, D., 2002, Accretion Power in Astrophysics, Cambridge University Press, Cambridge 
Fukue, J., 1983, PASJ, 35, 355 

Fukue, J., 1987, PASJ, 39, 309 

Fukue, J.& Yokoyama, T, Astronomical Society of Japan, Pubhcations (ISSN 0004-6264), vol. 40, no. 1, 1988, pp. 15-24 

Fukue, J., 2004, PASJ, 56, 681 

Fukue, J., 2004, PASJ, 56, 959 

Fukumura, K., Kazanas, D., 2007, ApJ, 669, 85 

Fukumura, K., & Tsuruta, S., 2004, ApJ, 611, 964 

Fukumura, K., Takahashi, M. & Tsuruta, S., 2007, ApJ, 657, 415 

Gammie, C. F, & Popham, R., 1998, ApJ, 498, 313 

Garfalo, D., 2009, ApJ, 699, 400 

Garfalo, D., 2009, ApJ, 699, L52 

Garfalo, D., Evans, D. A., & Sambuma, R. M., 2010, MNRAS, 406, 975 
Goswami, S., Khan, S. N., Ray, A. K., Das, T. K., 2007, MNRAS, 378, 1407 
Hawley J. R, & Krolik, J. H., 2001, APJ, 548, 348 

Hioki, M., & Maeda, K., 2009, Phys. Rev. D. vol. 80, Issue 2, id. 024042 Hioki, M., & Maeda, K., 2009, Phys. Rev. D. vol. 80, Issue 2, id. 
024042 

Ho, L. C. 1999, in Observational Evidence For Black Holes in the Universe, ed. S. K. Chakrabarti (Dordrecht: Kluwer), 153 



30 Das et. al. 



Huang, L., Cai, M., Shen, Z., & Yuan, R, 2007, MNRAS, 379, 833 
Hugoniot, H., 1887, Journal de I'Ecole Polytechnique 57, 3 
Hugoniot, H., 1887, Journal de I'Ecole Polytechnique 58, 1 
Igumenshchev, 1. V. & Abramowicz, M. A. 1999, MNRAS, 303, 309 
Igumenshchev, I. V., & Beloborodov, A. M., 1997, MNRAS, 284, 767 
lUarionov, A.R, & Sunyaev, R. A. 1975a, A & A, 39, 205 
niarionov, A. R 1988, Soviet Astron., 31, 618 

Janiuk, A., Sznajder, M., Moscibrodzka, M., & Proga, D., 2009, ApJ, 705, 1503 

Jordan, D. W., Smith, P., 1999, Nonlinear Ordinary Differential Equations, Oxford University Press, Oxford 
Kafatos, M., Yang, R. X., 1994, MNRAS, 268, 925 

Karas, V., Vokrouhlicky, D. & Polnarev, A. G., 1992, MNRAS (ISSN 0035-871 1), vol. 259, no. 3,pp. 569 

Kato, S., 1978, MNRAS, 185, 629 

Kato, S., Honma, R, Matsumoto, R., 1988, MNRAS, 231, 37 

Kato, S., Pukue, J., & Mineshige, S., 1998, Black Hole Accretion Disc, Kyoto University Press. 
Kato, Y, Miyoshi, M., Takahashi, R., Negoro, H., & Matsumoto, R., 2010, MNRAS, 403, L74 
Landau, L. D., Lifshitz, E. M., 1987, 'Pluid Mechanics', Butterworth-Heinemann, Oxford 
Landau, L. D., Lifshitz, E. M., 1994, 'Statistical Mechanics', Oxford: Pergamon 

Landau, L.D. & Lifshitz, E.M., 1981, Physical Kinetics, in Course of Theoretical Physics, Volume 10, Pergamon Press 

Lanzafame, G., 2008, ASPC, 385, 115 

Lemos, J. P S., & Letelier, P S., Phys. Rev. D. 49, 5135 

Liang, E. P T., 1984, Space. Sci. Rev. 38, 353 

Liang, E. P T., Thomson, K. A., 1980, ApJ, 240, 271 

Lu, J. R 1985, A & A, 148, 176 

Lu, J. R 1986, Gen. Rel. Grav. 18, 45L 

Lu, J. R, Yu, K. N., &. Young, E. C. M., 1995, A & A, 304, 662 

Lu, J. R, Yu, K. N., Yuan, R, Young, E. C. M., 1997, A & A, 321, 665 

Lu, J. R, Yu, K. N., Yuan, R, Young, E. C. M., 1997, Astrophysical Letters and Communications, 35, 389 

Lu, J. R, & Yuan, R, 1998, MNRAS, 295, 66 

Lu, J. R, & Gu, W. M., 2004, Chin. Phys. Lett., 21, 2551 

Luminet, J.-R, A & A, vol. 75, no. 1-2, May 1979, pp. 228 

Manmoto, T. 2000, ApJ, 534, 734 

Martinez-Sansigre, A., & Rawlings, S., 2011, MNRAS, 414, 1937 

Matsumoto, R., Kato, S., Pukue, J., Okazaki, A. T., 1984, PASJ, 36, 71 

McClintock, J. E., Narayan, R., Davis, S. W., Gou, L., Kulkami, A., Orosz, J. A., Penna, R. R, Remillard, R. A., & Steiner, J. P., 2011, 

Classical and Quantum Gravity, Voltmie 28, Issue 11, pp. 114009 
McNamara, B. R., Rohanizadegan, M., & Nulsen, P. E. J., 2011, ApJ, Volume 727, Issue 1, article id. 39 
Meliani, Z., Sauty, C, Tsinganos., K., & Vlahakis, N., 2004, A & A, Volume 425, pp. 773 
Migone, A. & McKinney, J.C., 2007, MNRAS, Volume 378, pp. 1118 

Miller, J. M., Reynolds, C. S., Rabian, A. C, Miniutti, G., & Gallo, L. C, 2009, ApJ, 697, 900 
Moscibrodzka, M., Das, T. K., & Czerny, B. 2006, MNRAS, 370, 219 
Moncrief, V., 1980, ApJ, 235, 1038 
Mondal, S., 2010, ApJ 708, 1442 

Mondal, S. & Basu, R, 2011, Class. Quant. Grav., 28, 235004 
Muchotrzeb, B., & Paczynski, B., 1982, Acta Actron. 32, 1 
Muchotrzeb, B., 1983, Acta Astron. 33, 79 
Muchotrzeb-Czerny, B., 1986, Acta Astronomica, 36, 1 
Mukhopadhyay, B., 1999, Indian J. Phys., Part B, Vol. 73B, No. 6, p. 917 
Mukhopadhyay, B. & Dutta, P., 2012, New Astronomy, 17, 51 

Nag, S., Acharya, S., Ray, A.K. & Das, Tapas K., 2012, New Astronomy, Volume 17, pp. 285 
Nagakura, H., Yamada, S., 2008, ApJ, 689, 391 
Nagakura, H., Yamada, S., 2009, ApJ, 696,2026 
Nakayama, K., Pukue, J., 1989, PASJ, 41, 271 
Nakayama, K., 1996, MNRAS, 281, 226 

Narayan, R., Mahadevan, R., & Quataert, E., 1998, in Theory of Black Hole Accretion Disks, ed. M. A. Abramowicz, G. Bjomsson, & J. 

E. Pringle (Cambridge; Cambridge Univ. Press), 343 
Nixon, C. J., Cossins, P J., King, A. R., & Pringle, J. E., 201 1, MNRAS 412, 1591 

Novikov, I., & Thome, K. S. 1973, in Black Holes,eds. c. De Witt and B. De Witt(Gordon and Breach, New York). 



Accretion close to event horizon 31 



Okuda, T., Teresi, V., Toscano, E., & Molteni, D., 2004, Publications of tlie Astronomical Society of Japan, Vol.56, No.3, pp. 547-552 

Okuda, T., Teresi, V., & Molteni, D., 2007, MNRAS, 377, 1431 

Paczynski, B., 1987, Nature, 327, 303 

Paczynski, B., Wiita R J., 1980, A&A, 88, 23 

Pariev, V. I., 1996, MNRAS, 283, 1264 

Parker, E. N., 1965, Space Science Reviews, vol. 4, p.666 

Peitz, J., Appl, S., 1997, MNRAS, 286, 681 

Popham, R., & Gammie, C. P., 1998, ApJ, 504, 419 

Pringle, J. E., 1981, ARA&A, 19, 137 

Proga, D., & Begelman, M. C, 2003, ApJ, 582, 69 

Rankine, W. J. M., 1870, Philosophical Transactions of the Royal Society of London 160: 277 

Ray, A. K., 2003a, MNRAS, 344, 83 
Ray, A. K., 2003b, MNRAS, 344, 1085 

Ray, A. K., Bhattacharjee, J. K., 2002, Phys. Rev. E, 66, 066303 

Ray, A. K., Bhattacharjee, J. K., 2007a, Classical and Quantum Gravity, 24, 1479 
Ray, A. K., Bhattacharjee, J. K., 2007b, Physics Letters A, 371, 241 

Reynolds, C. S, Brenneman, L. W., Lohflnk, A. M., Trippe, M. L., MiUer, J. M., Reis, R. C, Nowak, M. A., & Fabian, A. C., 2011, 'Probing 

Relativistic Astrophysics Around SMBHs: The Suzaku AGN Spin Survey', 
Riffert, H., & Herold, H., 1995, ApJ, 450, 508 

Ryu, D. & Chattopadhyay, L, 2006, Astrophys. J. Suppl. Series, 166(1), 410 

Salas, M. D., 2006, Shock Waves, 16, 477 

Sadowski, A., 2009, ApJS 183, 171 

Shakura, N. L, Sunyaev, R. A., 1973, A&A, 24, 337 

Sponholz, H., & Molteni, D., 1994, MNRAS, 271, 233 

Straub, O., Vincent, F. H., Abramowicz, M. A., Gourgoulhon, E., & Paumard, T., 2012, 'Modelling the black hole silhouette in Sgr A* with 

ion tori', arXiv: 12 03 . 2 618vl [astro-ph . GA] 
Strogatz, S., 2001, 'Nonlinear Dynamics And Chaos: With Applications To Physics, Biology, Chemistry, And Engineering', Westview 

Press; 1"* edition. 

Takahashi, M., Rilett, D., Fukumura, K., & Tsuruta, S., 2002, ApJ, 572, 950 
Takahashi, R., 2004, ApJ, 611, 996 

Takahashi, M., Goto, J., Fukumura, K., Rilett, D., & Tsuruta, S., 2006, ApJ, 645, 1408 
Takahashi, R., 2007, MNRAS, 382, 567 

Takahashi, M., Rilett, D., Fukumura, K., & Tsuruta, S., 1992, ApJ, 572, 950 
Tchekhovskoy, A., Narayan, R., & McKinney, J. C, 2010, ApJ, 71 1, 50 

Tchekhovskoy, A., & McKinney, J. C, 2012, MNRAS, Online early. DOL 10.1111/j.l745-3933.2012.01256.x 

Thorne, K. S., 1974, ApJ, 191, 507 

Toth, G., Keppens, R., & Botchev, M. A. 1998, A&A, 332. 1159 

Vincent, F. H., Gourgoulhon, E., Straub, O., Abramowicz, M., Novak, J., Paumard, T., Perrin, G., 2012, arXivl209.0251 

Wiita, R J., 2001, eprint (arXiv:astro-ph/0103020) 
Yang, R. X., Kafatos, M., 1995, A&A, 295, 238 

Yuan, F, Dong, S., & Lu, J. F, 1996, Astrophysics and Space Science, 246, 197 

Zakharov, A. F, Paolis, F. De., Ingrosso, G., & Nucita, A. A., 2012, New Astronomy Review, Volume 56, Issue 2, p. 64-73 
Ziolkowski, J., 2010, Memorie della Societ Astronomica Italiana, v.81, p.294 



